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INTRODUCTION 


1.1  STATEMENT  OF  THE  PROBLEM 

The  ability  to  monitor  and  predict  the  intensity"*-  of 
tropical  cyclones  is  important  because  of  the  death  and 
destruction  caused  by  the  high  wind  speeds,  storm  surges, 
and  heavy  precipitation  which  characterize  these  cyclones 
(Anthes,  1982) .  Tropical  cyclones  occur  over  the  warm 
tropical  oceans  where  they  usually  affect  coastal  habita¬ 
tions  and  ships  at  sea.  Because  of  the  sparse  synoptic  re¬ 
porting  network  in  these  areas,  an  important  source  of 
tropical  cyclone  intensity  data  is  reconnaissance  by  air¬ 
craft.  However,  in  recent  years  as  aircraft  reconnaissance 
has  been  reduced  for  economic  reasons,  an  effort  has  been 
made  to  fill  the  resulting  data  gaps  by  using  meteorologi¬ 
cal  satellite  information  (Gentry  et  al . ,  1980). 

Among  all  meteorological  satellites,  the  geostationary 
satellites  are  particularly  suited  for  tropical  cyclone 
data  collection.  Because  they  image  the  same  disk  of  the 
Earth  every  half-hour  they  acquire  data  more  often  than  is 
feasible  using  aircraft  reconnaissance.  The  usable  data 
collected  are  the  infrared  ( IR)  radiance  count  values  which 
can  be  converted  to  equivalent  blackbody  cloud  top  tempera¬ 
tures  (Gentry  et  al.,  1980). 

"*"Defined  as  the  maximum  sustained  surface  wind  speed. 
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In  the  past,  the  satellite  information  has  been  inter¬ 
preted  in  terms  of  the  spiral-shaped  cloud  patterns  of 
tropical  cyclones  (Arnold,  1977)  .  Dvorak  (19  75)  empirically 
related  the  patterns  to  tropical  cyclone  intensity  by  using 
visible  and  IR  satellite  images.  However,  with  the  develop¬ 
ment  of  interactive  computer  techniques  for  manipulating 
satellite  data,  it  has  become  evident  that  the  relationship 
of  the  cloud  patterns  in  the  digital  IR  data  to  tropical 
cyclone  intensity  can  be  better  quantified. 

1.2  OBJECTIVE 

The  specific  purpose  of  this  study  is  to  determine  a 
quantitative  relationship  between  satellite  measurements  of 
the  IR  radiance  field  and  the  operationally  reported  inten¬ 
sity  of  tropical  cyclones. 

1 . 3  APPROACH 

The  approach  to  the  problem  consists  of  three  steps: 

i)  A  detailed  survey  of  past  studies  to  gain  an  insight 
into  the  relationship  of  observed  spiral-shaped  cloud 
patterns  to  tropical  cyclone  intensity. 

ii)  The  development  of  the  methodology  to  apply  the  past 
research  findings  to  extract  the  relevant  parameters  from 
the  digital  IR  radiance  field. 

iii)  The  development  and  testing  of  multiple  linear  re¬ 
gression  equations  for  monitoring  and  predicting  tropical 
cyclone  intensity  using  the  extracted  parameters. 
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2.  A  SURVEY  OF  PAST  STUDIES  OF  THE  RELATIONSHIP  OF 
CLOUD  PATTERNS  TO  TROPICAL  CYCLONE  INTENSITY 

2.1  OBSERVED  SPIRAL  PATTERNS 

In  studies  of  radar  films  of  tropical  cloud  features 
Senn  and  Hiser  (1959)  found  that  hurricanes  are  character¬ 
ized  by  asymmetric  rain  bands  that  generally  conform  to 
logarithmic  spirals.  The  equation  of  a  logarithmic  spiral 
is 


£n  r  =  -  a  0  (1) 

where  r  is  the  radius,  a  is  the  inflow  or  crossing  angle 
that  the  spiral  makes  with  concentric  circles  about  the 
spiral  center,  and  0  is  the  azimuth.  They  found  that  these 
spiral  rain  bands  usually  extended  100  to  200  km  from  the 
cyclone  center,  and  that  there  were  often  several  spiral 
bands  associated  with  each  cyclone.  The  spiral  cloud  bands 
were  composed  of  numerous  large  cellular  convective  clouds 
that  propagated  along  the  band  towards  the  spiral  center. 

Lahiri  (1981)  used  a  simple  model  of  a  tropical  cyclone 
in  which  the  low-level  streamlines  were  logarithmic  spirals 
and  showed  that  the  rate  of  generation  of  latent  heat  in  the 
model  was  proportional  to  a.  Lahiri  observed  that  many 
satellite  images  of  tropical  cyclones  exhibit  cloud  banding 
patterns  that  conform  to  logarithmic  spirals  and  assumed 
that  the  cloud  band  spirals  were  coincident  with  the 
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streamline  spirals.  Lahiri  then  showed  that  the  spiral 
cloud  bands  degenerated  to  circles  (a  approached  zero)  as 
the  cyclones  matured. 

The  relationship  between  latent  heating  and  cloud  pat¬ 
terns  is  important  because  the  latent  heat  released  within 
deep  cumulus  towers  is  the  tropical  cyclone's  primary  energy 
source  (Anthes,  1982) .  The  satellite  senses  the  IR  radi¬ 
ance  field  of  the  tropical  cyclone  and  its  environment.  The 
radiance  values  for  the  cloud  tops  can  be  converted  to 
equivalent  blackbody  cloud  top  temperatures  which  can  pro¬ 
vide  quantitative  information  on  the  availability  of  latent 
heat  within  the  cyclone  (Gentry  et  al.,  1980).  Therefore, 
the  latent  heat  released  in  a  cyclone  is  indicated  by  the 
amount  and  strength  of  the  convection,  which  are  related  to 
the  areal  distribution  of  the  cyclone's  satellite-derived 
IR  radiance  field. 

In  a  linear  stability  analysis  of  the  development  of 
spiral  cloud  bands  in  tropical  cyclones  Kurihara  (1976) 
found  that  the  preferred  spiral  scale  of  the  large  scale 
bands  was  wavenumber  two  with  a  band  width  of  200  km.  The 
flow  perturbation  fields  of  spiral  shape  were  modeled  by 
logarithmic  spirals  with  modifying  Bessel  functions. 
Willoughby  (1978)  showed  a  linear  model  of  hurricane  rain 
bands  in  cylindrical,  log-pressure  coordinates  that  gener¬ 
ated  a  physically  realistic  radial  geopotential  structure 
with  waves  of  wavenumber  two  and  higher . 
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Black  and  Anthes  (1971)  quantified  the  spiral  structure 
of  the  tropical  cyclone  outflow  layer  by  harmonic  analysis 
of  the  tangential  and  radial  components  of  the  satellite 
derived  flow  fields.  Their  analysis  of  five  tropical  cy¬ 
clones  showed  the  most  significant  number  of  spiral  bands 
to  be  less  than  four  and  that  the  horizontal  scale  of  the 
bands  increased  outward  from  the  storm  center.  In  a  more 
recent  study,  Rodgers  and  Gentry  (1983)  traced  the  upper 
and  lower  tropospheric  wind  fields  around  three  tropical 
cyclones  from  satellite-derived  cloud  drift  winds.  They 
found  that  the  flow  at  both  levels  revealed  a  large  scale 
spiral  structure.  They  used  the  fields  to  compute  the 
local  changes  in  the  cyclone's  net  relative  angular  momentum 
and  found  that  it  appeared  to  be  related  to  the  trend  in 
storm  intensity. 

2.2  SATELLITE-DERIVED  ESTIMATES  AND  PREDICTIONS  OF 

TROPICAL  CYCLONE  INTENSITY 

Techniques  for  estimating  and  predicting  the  intensity 
of  tropical  cyclones  from  satellite  IR  or  visible  pictures 
focus  on  the  identification  of  various  tropical  cloud 
features,  such  as  the  shape  and  areal  extent  of  the  cirrus 
overcast,  the  appearance  of  the  eye,  and  the  amount  of  cloud 
banding  present  (Fett,  1966;  Fritz  et  al.,  1966;  Hubert  et 
al.,  1969;  Dvorak,  1973,  1975,  and  1979).  The  widely  used 
Dvorak  (1975)  method  is  recognized  as  the  most  accurate 
subjective  technique  for  estimating  the  current  and  future 
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intensity  of  tropical  cyclones  using  satellite  imagery 
(Shewchuk  and  Weir,  1980).  The  basic  premise  of  Dvorak's 
technique  is  that  by  associating  the  organization  of  the 
various  cloud  features  of  tropical  cyclones  with  their 
intensities,  estimates  of  the  current  and  future  intensi¬ 
ties  for  similarly  organized  cyclones  can  be  obtained. 

Figure  1  shows  the  characteristic,  satellite-observed 
cloud  features  that  are  related  to  tropical  cyclone  inten¬ 
sity  by  Dvorak's  method.  The  central  cloud  features  (CF) 
refer  to  the  extent  and  shape  of  the  cloud  elements  that 
appear  within  the  curve  of  the  comma  band  and  can  either 
surround  or  cover  the  storm  center.  Intensification  is 
indicated  by  the  CF  becoming  colder,  or  by  the  eye  becoming 
either  warmer,  smaller,  or  more  circular.  The  cyclone  s 
outer  banding  features  (BF)  in  Figure  1  refer  to  the  outer 
comma  shaped  cloud  band  that  curves  evenly  around  the 
central  features.  Intensification  is  indicated  by  the  BF 
becoming  more  developed. 

In  the  application  of  the  technique  the  CF  and  BF  are 
evaluated  and  their  results  combined  to  provide  a  "T-number" 
description  of  the  intensity  of  a  cyclone.  The  assigned  T- 
number  is  plotted  on  a  development  diagram  (see  Fig.  2)  and 
the  cyclone  is  expected  to  intensify  or  weaken  according  to 
one  of  the  three  curves  plotted.  Rapidly  developing  or 
rapidly  weakening  cyclones  (±1.5  T— numbers  per  day)  are 
plotted  on  the  RAPID  curve,  cyclones  exhibiting  normal 
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Figure  1.  A  schematic  of  a  typical  IR  satellite  image 
showing  the  different  cloud  features  which  are  related 
to  tropical  cyclone  intensity  by  Dvorak's  (1979)  method. 
The  contours  are  isolines  of  equal  gray  shades.  The 
contour  values  increase  towards  the  cyclone  center  and 
are  inversely  proportional  to  the  temperatures  of  the 
cloud  tops.  The  circular  area  in  the  center  of  the  cy¬ 
clone  is  the  warm  eye.  The  area  inside  of  the  dashed 
lines,  the  cloud  features  encircling  the  eye,  are  the 
central  features  (CF) .  The  cloud  bands  outside  of  the 
dashed  lines  are  the  banding  features  (BF) . 
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Figure  2.  Intensity  change  curves  of  the  Dvorak  (1975) 
model.  The  T  (for  tropical)  numbers  are  defined  in 
terms  of  the  cloud  features  of  the  storm  and  relate  to 
the  storm's  intensity.  The  curves  depict  tropical  cyclone 
development  and  weakening  as  occurring  along  one  of 
the  three  curves  of  T-number  vs.  time.  A  typical  24  h 
change  is  one  T-number  (Dvorak,  1975) .  The  relationship 
between  the  T-numbers  and  the  maximum  sustained  surface 
wind  speed  is  shown  to  the  left  in  the  figure. 
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developing  or  weakening  (±1.0  T-numbers  per  day)  are  plot¬ 
ted  on  the  TYPICAL  curve,  and  slowly  developing  or  slowly 
weakening  cyclones  (±0.5  T-numbers  per  day)  are  plotted  on 
the  SLOW  curve.  A  24  h  forecast  is  made  by  incrementing 
the  abscissa  one  day  along  the  cyclone's  characteristic 
curve.  The  C.I.  number  (current  intensity)  on  the  WEAKENING 
portion  of  the  figure  is  a  modified  T-number  to  account  for 
factors  which  are  not  directly  related  to  the  BF  or  CF  (see 
Dvorak,  1975) .  The  relationship  between  T-number  and  maxi¬ 
mum  surface  wind  speed  can  be  seen  in  the  left  of  the  fig¬ 
ure.  Note  that  the  larger  T-numbers  encompass  a  wide  range 
of  wind  speeds. 

An  objective  method  of  predicting  tropical  cyclone  in¬ 
tensity  by  using  computer  processed  IR  satellite  data  has 
been  developed  by  Gentry  et  al.  (1980)  .  Radially  averaged 
equivalent  blackbody  cloud  top  temperatures  (derived  from 
the  IR  radiance  field)  about  the  cyclone  center  and  their 
standard  deviations  were  used  to  make  24  h  predictions  of 
tropical  cyclone  intensity  via  linear  regression  equations. 
They  found  that  the  relationship  of  the  IR  data  to  future 
intensity  was  nonlinear  over  time.  Consequently,  their 
best  results  were  reported  when  the  data  were  stratified 
based  on  the  current  intensity,  which  was  estimated  by  using 
the  Dvorak  (1975)  technique.  The  tropical  cyclones  were 
stratified  into  two  groups:  current  intensity  less  than  or 
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equal  to  65  kt  and  current  intensity  greater  than  65  kt. 
Gentry  et  al.  (1980)  concluded  that  the  areal  distribution 
of  the  IR  radiance  field  serves  as  an  index  of  the  organi¬ 
zation  of  the  cyclone's  convective  activity  and  that  the 
relationship  can  be  successfully  exploited  to  predict 
tropical  cyclone  intensity. 

Swadley  (1983)  used  a  spiral  linearization  technique 
( SLT)  on  an  interactive  computer  to  transform  IR  satellite 
images  of  the  cloud  bands  of  tropical  cyclones  into  loga¬ 
rithmic  spiral  coordinates.  Using  a  limited  number  of 
satellite  images,  Swadley  developed  linear  regression  equa¬ 
tions  for  estimating  tropical  cyclone  intensity  based  on  the 
spiral  inflow  angle  (a)  and  three  measurements  of  the  IR 
radiance  field.  He  found  that  the  equations  could  not  ex¬ 
plain  the  observed  variability  of  intensity  between  the 
cyclones  and  concluded  that  more  information  from  the  digi¬ 
tal  IR  satellite  data  was  necessary  in  order  to  develop  a 
better  relationship  with  tropical  cyclone  intensity. 

Using  a  large  number  of  satellite  images  processed  by 
the  SLT,  Peterson  et  al .  (1983)  evaluated  the  statistical 

relationship  between  a  and  tropical  cyclone  intensity. 

They  showed  that  a  behaved  randomly  with  time,  contrary  to 
Lahiri's  (1981)  findings,  and  concluded  that  it  was  impos¬ 
sible  to  characterize  the  intensity  of  a  tropical  cyclone 
using  a  single  parameter. 
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In  summary,  past  research  has  shown  that  tropical  cy¬ 
clone  cloud  bands  are  characterized  by  a  low  wavenumber 
spiral  pattern.  The  spiral-shaped  cloud  bands  are  observed 
in  IR  satellite  data  and  have  been  empirically  related  to 
tropical  cyclone  intensity. 

Recently,  satellite  data  manipulation  techniques  on 
interactive  computers  have  shown  promise  for  better  quanti— 
fying  the  relationship  of  the  cloud  bands  to  tropical  cy¬ 
clone  intensity.  The  method  that  shows  the  most  promise  is 
spiral  analysis  coupled  with  multiple  linear  regression. 
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3.  DEVELOPMENT  OF  METHODOLOGY 


In  the  present  chapter  linear  regression  techniques  are 
used  to  establish  a  relationship  between  the  spiral-shaped 
cloud  bands  and  the  intensity  of  tropical  cyclones.  A  major 
step  in  the  process  is  the  development  of  a  spiral  analysis 
technique  to  quantify  the  spiral  patterns  of  the  cloud 
bands.  In  this  method  the  spiral  characteristics  are  ex¬ 
tracted  from  the  satellite-derived  digital  IR  data  using 
multiple  Fourier  analyses.  The  Fourier  analyses  are  per¬ 
formed  over  a  spiral  grid  which  is  unique  to  each  image  of 
a  tropical  cyclone. 

This  particular  technique  was  chosen  because  it  incor¬ 
porates  the  large  scale  spiral  shape  of  the  cloud  bands  into 
the  Fourier  analyses.  The  technique  differs  from  previous 
techniques  (e.g. ,  Swadley ,  1983)  in  that  it  is  performed  in 
geographic  coordinates  instead  of  local  polar  coordinates, 
and  the  data  are  analyzed  in  place  instead  of  being  trans¬ 
formed  into  spiral  coordinates.  Therefore,  there  is  no 
distortion  due  to  the  curvature  of  the  earth  or  coordinate 
transforms . 

3.1  SPIRAL  ANALYSIS  TECHNIQUE 

Each  digital  IR  satellite  image  processed  in  the  pres¬ 
ent  study  was  analyzed  on  the  Satellite— data  Processing  and 
Display  System  (SPADS)  Eclipse  S-250  minicomputer  at  the 
Naval  Environmental  Prediction  Research  Facility  (NEPRF) . 
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The  unique  capabilities  of  the  SPADS  for  interactive  digi 
tal  image  processing  are  described  by  Schramm  et  al . 

(1982)  . 

The  relevant  subroutines  of  the  computer  program  (TEDS) 
written  for  this  study  to  implement  the  spiral  analysis 
technique  on  the  SPADS  are  reproduced  in  Appendix  A.  TEDS 
utilizes  an  algorithm  adapted  from  Nicholson  (1983)  to  de¬ 
termine  the  spiral  analysis  grid  for  an  image  and  to  per¬ 
form  the  Fourier  analyses  of  the  IR  satellite  data  over  the 
spiral  grid. 

3.1.1  Determination  of  the  Spiral  Analysis  Grid 

Figure  3  is  an  example  of  a  typical  spiral  analysis 
grid.  The  grid  is  defined  as  ten  parallel  spirals  of  data 
points  which  have  a  common  center  and  inflow  angle  (a) .  The 
spirals  are  separated  from  each  other  by  18°  and  are  or¬ 
thogonal  to  the  cloud  bands.  The  inflow  angle  and  spiral 
center  of  the  spiral  analysis  grid  are  determined  by  an 
average  of  the  spiral  fits  to  all  of  the  cloud  bands  proc¬ 
essed  on  a  particular  image.  The  processing  consists  of 
calculating  the  best  least  squares  fit  spiral  corresponding 
to  each  cloud  band  and  then  averaging  them.  The  detailed 
procedure  is  as  follows:  after  displaying  the  image  on 
the  video  monitor,  the  image  coordinates  (lines,  elements) 
of  the  cloud  band  are  entered  into  TEDS  via  the  intercon¬ 
nected  graph  pen/tablet.  The  pen  is  drawn  along  the  center 
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Figure  3.  The  dots  represent  a  typical  spiral  analysis 
grid  used  in  this  study.  Each  dot  represents  a  data 
point.  The  grid  consists  of  ten  parallel  spiral  rows  of 
data  points  which  have  a  common  center  (not  shown)  and 
are  orthogonal  to  the  cloud  bands.  Each  spiral  has  the 
same  inflow  angle  a  and  is  separated  from  its  neighbors 
by  18°.  How  the  grid  is  determined  and  oriented  on  an 
image  is  described  later  in  the  text. 
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of  the  cloud  band,  parallel  to  the  edge,  starting  at  its 
outermost  point  and  working  inwards  toward  the  center 
( see  Fig .  4 ) . 

The  entered  coordinates  are  earth-located  (navigated) 
and  a  spherical— logarithmic  spiral  is  least  squares  fitted 
to  them  following  Nicholson  (1982) .  This  process  is 
thoroughly  described  in  Appendix  B.  Figure  5  is  a  mercator 
plot  of  a  spherical-log  spiral  and  its  representation  in 
spiral  coordinates.  The  equation  of  a  spherical-log  spiral 
is 

y  -tana(v-vQ) 

y  =  2  arctan  ( tan  (~2~)  e  )  ,  (2) 

where  y  is  the  great  circle  distance,  and  v  the  angular 

separation,  between  points  along  the  spiral  and  the  spiral 

center;  y  and  v  refer  to  the  same  quantities  at  the  outer- 
o  <o 

most  point  along  the  spiral  and  a  is  the  inflow  angle  of 
the  spiral. 

Since  there  are  three  unknowns  (a,  yQ,  and  vq)  and  the 
least  squares  formulation  yields  only  two  normal  equations, 
an  iterative  computational  scheme  is  used  to  calculate  one 
variable.  The  normal  equations  are  most  easily  solved  for 
a  and  v  ,  therefore,  yQ  is  determined  using  a  highly  effi¬ 
cient  iterative  computational  scheme  following  Nicholson 
(1982) .  The  iterative  scheme  is  used  to  converge  on  the 
best  fit  spiral  center,  which  uniquely  determines  yQ. 
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Figure  4.  The  heavy  line  shown  with  starting  and  ending 
points  is  the  tracing  of  the  graph  pen  along  the  center 
of  the  schematically  represented  cloud  band.  In  this 
example,  the  cloud  band  can  be  traced  no  farther  towards 
the  cyclone  center  (eye) . 
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Figure  5.  The  inset  is  a  mercator  plot  of  the  data  points 
entered  along  the  heavy  line  shown  in  Fig.  4.  An  X  marks 
the  best  fit  spiral  center  and,  thus,  defines  yQ  and  vQ. 
The  equation  of  a  spherical-log  spiral  is 
y  =  2  arctan  [tan  (\iq/2)  exp  (-  tan  a  ( v-vQ)  )  ]  •  An  arc  of  a 
circle  about  the  spiral  center  is  shown,  defining  a,  which 
is  41°.  The  bottom  graph  are  the  data  points  shown  in 
the  inset  plotted  in  spherical- log  spiral  coordinates. 

The  representation  of  a  in  the  spiral  coordinates  is  also 
shown . 
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The  upper  left  and  lower  right  corners  of  the  area  in 
which  to  search  for  the  spiral  center  are  entered  from  the 
graph  tablet  (see  Fig.  6) .  The  search  area  is  divided  into 
a  six-by-six  grid  and  each  interior  grid  point  used  as  a 
sample  spiral  center.  The  following  solutions  to  the  normal 
equations  (Appendix  B)  are  then  solved  for  vq  and  a: 


cot  a  = 


N  N  N 

£  V.  1  Sin  T.  -  N  £  V  .  in  T. 

1=1  1=1 _ 1=1 _ 

N  ~  N  2 

N  l  (£n  T. )  -  (  l  Hn  T. ) 

i=l  1  i=l  1 


(3) 


where 


N 


N 


N 


l  V  l  Un  T  )  -  l  £n 

i=l  1  i=l  1  i=l 


N 

I 


v  .  £n  T  . 
1  i=l  1  1 


N  2  N  ? 

N  l  (£n  T.)  -  (  l  £n  T  ) 

i=l  1  i=l  1 


(4) 


T  . 
1 


.  ,Vi, 

tan  (-2-) 

tan  (-^-) 


The  iterative  scheme  continues  as  the  sample  spiral 
center  with  the  smallest  root  mean  square  error  is  focused 
on  as  the  new  grid  center  and  the  grid  size  divided  by  six. 
Figure  7  is  a  schematic  of  the  process,  showing  three  of  the 
nine  iterations  used  to  bring  the  grid  spacing  to  image 
resolution.  This  minimizes  the  RMS  error  while  determining 
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Figure  6.  The  upper  left  and  lower  right  coordinates  of  the 
area  in  which  to  search  for  the  spiral  center  of  the  entered 
heavy  line  are  entered  via  the  graph  pen/tablet.  The 
drawn  box  is  a  spherical  quadrilateral  whose  coordinates 
are  calculated  and  displayed  by  TEDS.  The  size  of  the 
box  is  arbitrary. 


THIS  POINT  ENTERED 
FIRST  USING 
THE  GRAPH 
PEN/TABLET. 


THIS  POINT 
ENTERED  SECOND 
USING  THE 
GRAPH  PEN/ 
TABLET. 
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Figure  7.  A  mercator  representation  of  the  box  described  in 
Fig.  6.  The  triangles  at  the  corners  represent  the  user 
entered  points.  The  box  is  divided  into  a  6x6  grid  and 
the  25  interior  points  (circled)  are  used  as  sample  spiral 
centers  and  their  RMS  error  calculated.  Assuming  the 
point  II, J4  has  the  smallest  RMS  error,  it  is  focused  on 
as  the  new  grid  center  and  the  grid  size  is  divided 
further  by  six.  The  process  is  repeated  on  the  new  grid 
(solid  dots),  and  assuming  the  point  K5,L4  has  the  smallest 
RMS  error,  it  is  focused  on  as  the  new  grid  center  and 
the  grid  size  divided  by  six  again.  The  process  is  repeated 
nine  times  to  bring  the  grid  spacing  to  image  resolution. 
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a,  iaQ/  vq,  and  the  latitude  and  longitude  of  the  best  fit 
spiral  center. 

The  five  parameters  calculated  from  the  input  line 
along  the  cloud  band  are  then  used  to  reconstruct  the  best 
fit  spherical-log  spiral  on  the  video  monitor.  Using  the 
plotted  spiral  fit  and  RMS  error  statistic  as  data  quality 
checks,  the  operator  then  decides  whether  or  not  to  save 
the  calculated  data  for  use  in  the  determination  of  the 
spiral  analysis  grid.  During  processing,  the  inflow  angles 
and  spiral  centers  saved  by  the  operator  are  averaged.  As 
noted  earlier,  the  averaged  quantities  determine  the  spiral 
characteristics  of  the  analysis  grid. 

The  final  step  in  the  process  of  determining  the  spiral 
analysis  grid  is  orienting  it  properly  on  the  image  so  that 
the  desired  cloud  pattern  is  Fourier  analyzed.  The  orien¬ 
tation  is  accomplished  by  defining  a  delimiting  spiral  with 
the  graph  pen.  As  shown  in  Figure  8,  the  delimiting  spiral 
has  the  averaged  center  and  inflow  angle  described  above,  how¬ 
ever,  it  is  rotated  with  respect  to  the  cloud  pattern  to  be 
analyzed:  and  vq  are  selected  such  that  the  cloud  pat¬ 

tern  of  interest  is  located  between  the  spiral  center  and 
the  delimiting  spiral. 

An  important  point  is  that  in  spiral  coordinates  the 
delimiting  spiral  is  orthogonal  to  the  ten  spirals  which 
comprise  the  analysis  grid.  The  delimiting  spiral  only 
defines  the  orientation  of  the  grid  on  the  image. 
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Figure  8.  The  heavy  line  is  the  delimiting  spiral  of  the 
tropical  cyclone  and  defines  the  orientation  of  the  spiral 
analysis  grid.  The  latitude  and  longitude  of  its  center, 
and  its  inflow  angle  a,  are  those  of  the  average  of  all 
of  the  cloud  bands  processed.  y  and  v  are  selected  such 
that  the  cloud  band  to  be  analyzed  is  l8cated  between  the 
spiral  center  and  the  delimiting  spiral.  v  is  measured 
counterclockwise  from  a  latitude  circle.  ° 
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3.1.2  Spiral-Fourier  Analysis  of  the  IR  Radiance  Field 

In  order  to  quantify  the  spiral  patterns  of  the  cloud 
bands  multiple  spiral-Fourier  analyses  (Fourier  analyses 
along  spherical-log  spirals)  are  used.  The  analyses  are 
performed  over  the  spiral  grid,  along  the  ten  spirals 
orthogonal  to  the  delimiting  spiral  (see  Fig.  9) .  There¬ 
fore,  the  Fourier  analyses  decompose  the  IR  radiance  field 
into  spiral  harmonics  parallel  to  the  delimiting  spiral 
(see  Fig.  10) . 

The  spiral-Fourier  analyses  proceed  as  follows:  along 
each  of  the  orthogonal  spirals  of  the  grid  the  geographic 
locations  of  forty  equally  spaced  data  points  are  calcu¬ 
lated.  In  spiral  coordinates  distances  are  measured  as  a 
function  of  the  natural  logarithm  of  the  radius,  therefore, 
data  points  separated  by  equal  intervals  physically  become 
farther  apart  as  they  are  sampled  farther  from  the  spiral 
center.  As  the  data  are  sampled  closer  to  the  spiral  center 
(as  the  radius  goes  to  zero)  they  become  indistinguishable 
on  the  finite  resolution  image.  In  order  to  avoid  this 
difficulty,  data  are  not  sampled  along  an  orthogonal  spiral 
until  the  orthogonal  spiral  has  passed  beyond  the  inner 
portion  of  the  delimiting  spiral  (see  Figs.  9  and  10). 

After  the  locations  of  the  forty  data  points  along  the 
orthogonal  spiral  are  calculated,  their  IR  radiance  values 
are  retrieved  and  their  conventional  equal  spaced  discrete 
Fourier  transform  is  calculated  (Panofsky  and  Brier,  1968) . 
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DELIMITING  SPIRAL 
FROM  FIG.  8. 
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Figure  9.  The  spiral  analysis  grid  for  Fig.  8.  The 

numbered  dotted  lines  represent  the  ten  spiral  sections 
along  which  the  spiral-Fourier  analyses  of  the  IR  radi¬ 
ance  field  are  taken.  The  dotted  lines  are  orthogonal 
to  the  delimiting  spiral  and  are  composed  of  40  equal 
spaced  data  points  each.  Note  that  the  data  points  are 
not  sampled  until  after  each  orthogonal  spiral  has  passed 
over  the  delimiting  spiral.  This  is  necessary  so  that 
the  same  pixels  are  not  sampled  more  than  once. 
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. . —  HYPOTHETICAL  WAVE  FORM 

OF  WAVE  NUMBER  20. 


Figure  10 .  The  parallel  spirals  show  how  the  Fourier  decom¬ 
position  orthogonal  to  the  delimiting  spiral  resolves 
spiral  harmonics.  The  spirals  shown  are  for  wavenumber 
20,  which  is  represented  in  part  by  the  sine  wave.  The 
actual  analysis  grid  of  Fig.  9  extends  between  the  spiral 
sections  1  and  10. 


Preliminary  analyses  showed  that  in  some  cases  the 
variance  was  shifted  artificially  into  the  high  wavenumbers. 
This  was  due  to  distortion  (spikes)  in  the  data:  data 
sampled  containing  the  radiance  value  of  a  point  at  the 
surface  of  the  Earth  (relatively  warm)  instead  of  that  of  a 
cloud  top  (relatively  cold) .  Because  of  the  nearly  constant 
spectrum  with  respect  to  a  large  mean,  the  spikes  consti¬ 
tute  high  frequency  waves  with  large  amplitudes  that 
contribute  to  a  large  variance. 

Ridding  the  data  of  large  amplitude  distortion  is 
referred  to  as  prewhitening  since  the  ideal  is  to  bring 
the  spectrum  close  to  that  of  white  noise.  According  to 
Blackman  and  Tukey  (1958)  the  flattening  of  the  spectrum 
need  not  be  precise  since  we  need  only  to  make  the  change 
of  the  power  spectrum  with  frequency  relatively  small.  The 
methodology  developed  for  the  present  study  to  prewhiten 
the  data  is  to  identify  the  distortions  in  the  data  and  set 
them  equal  to  -1.5  standard  deviations  from  the  mean.  A 
lower  threshold  IR  value  is  calculated  for  each  of  the  ten 
spiral  sections  of  the  analysis  using  the  above  method, 
and  the  data  prewhitened,  prior  to  the  Fourier  analyses. 
Since  only  extremely  low  values  of  radiances  are  affected 
by  this  procedure,  approximately  95%  of  the  data  are  passed 
unaffected. 

2The  forty  data  points  calculated  for  each  orthogonal 
spiral  are  referred  to  as  spiral  sections  or,  if  there  is 
no  ambiguity,  simply  as  sections. 
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3.2  STATISTICAL  ANALYSIS 

Regression  equations  of  the  form 

y  =  a  +  fc>1x1  +  b2x2  +  ...  +  bpxp  (5) 

are  used  to  relate  spiral  parameters  and  past  intensity 
measurements  to  current  and  future  tropical  cyclone  inten¬ 
sity.  The  symbols  are: 

y  are  tropical  cyclone  intensity  estimates  at  0,  12, 
and  24  h  from  image  time, 

x  ,...,x  are  the  independent  variables  screened, 
i  P 

b  , . . . ,b  are  the  estimated  coefficients, 

1  P 

a  is  the  estimated  intercept,  and 

p  is  the  number  of  independent  variables  picked. 

The  hypothesis  is  that  the  intensity  of  tropical  cy¬ 
clones  is  related  to  their  spiral  banding  structure  as 
discussed  in  Chapter  2.  In  order  to  test  this  hypothesis, 
the  past  intensity  measurements  and  the  results  of  the 
spiral-Fourier  analyses  must  be  reduced  to  an  appropriate 
field  of  potential  predictors  for  the  stepwise  linear  re¬ 
gression  program. 

3.2.1  Selection  of  Potential  Predictors 

Dvorak  (1975)  has  shown  that  persistence  is  an  impor¬ 
tant  predictor  for  tropical  cyclone  intensity.  Therefore, 
tropical  cyclone  intensities  24  and  12  h  before  image  time. 


27 


and  their  difference,  are  used  as  possible  predictors  (see 
Table  1) .  For  the  intensity  forecast  algorithms  the  inten¬ 
sity  changes  over  the  past  24  and  12  h  using  the  estimated 
current  intensity  are  calculated  and  used  as  possible 
predictors. 

Swadley  (1983)  and  Lahiri  (1981)  showed  that  the  in¬ 
flow  angle  (a)  of  the  spiral  cloud  bands  was  also  an  impor¬ 
tant  intensity  predictor.  To  include  the  possibility  of 
that  relationship  in  this  study,  a  for  the  spiral  analysis 
grid  is  also  used  as  a  potential  predictor. 

Table  1  also  shows  that  additional  potential  predic¬ 
tors  are  considered.  The  position  of  the  cyclone  at  image 
time  is  used  to  include  possible  latitudinal  and  longi¬ 
tudinal  differences  in  intensity.  Diurnal  effects  are 
parameterized  using  the  image  time  and  seasonal  effects  are 
parameterized  using  the  Julian  date  of  the  image. 

The  large  amount  of  data  calculated  using  the  spiral- 
Pourier  analyses  is  reduced  by  limiting  the  contribution  to 
the  field  of  potential  predictors  to  representations  of  the 
spiral  signal  and  low  wavenumber  characteristics  of  the  IR 
radiance  field.  This  procedure  is  in  agreement  with  the 
results  of  Black  and  Anthes  (1971),  Kurihara  (1976),  and 
Willoughby  (1978)  who  showed  that  the  spiral  flow  fields  of 
tropical  cyclones  exhibit  large  signal  to  noise  ratios  for 
the  low  wavenumbers.  As  shown  in  Table  1,  the  potential 
predictors  chosen  to  represent  the  low  wavenumber 
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Table  1.  Potential  predictors 


REPRESENTATION 

DESCRIPTION 

LOW  WAVENUMBER 
CHARACTERISTICS 

Number  of  harmonics  necessary  to 
explain  85%  of  the  observed  variance 
for  sections  2,5,8,  and  the  mean  of 
all  sections. 

Number  of  harmonics  necessary  to  ex¬ 
plain  90%  of  the  observed  variance 
for  sections  2,5,8,  and  the  mean 
of  all  sections. 

Number  of  harmonics  necessary  to 
explain  95%  of  the  observed  variance 
for  sections  2,5,8,  and  the  mean  of 
all  sections. 

SPIRAL  SIGNAL 

Sum  of  the  variance  explained  by  the 
first  three  harmonics  for  sections 
2,5,8,  and  the  mean  of  all  sections. 

Number  of  the  harmonic  explaining  the 
maximum  variance  for  sections  2,5,8, 
and  the  mean  of  all  sections. 

Amplitude  of  the  harmonic  explaining 
the  maximum  variance  for  sections 
2,5,8,  and  the  mean  of  all  sections. 

VARIANCE 

Total  observed  variance  of  sections 
2,5,8,  and  the  mean  of  all  sections. 

Standard  deviation  of  sections  2,5, 

8,  and  the  mean  of  all  sections. 

Threshold  value  of  sections  2,5,8, 
and  the  mean  of  all  sections. 
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Table  1. (CONTINUED) 


REPRESENTATION 

DESCRIPTION 

PERSISTENCE 

Intensity  24  h  before  image  time. 

Intensity  12  h  before  image  time. 

Intensity  change  from  -24  h  to  -12  h 
before  image  time. 

Intensity  change  from  -24  h  to  image 
time  (using  predicted  current  inten¬ 
sity)  . 

Intensity  change  from  -12  h  to  image 
time  (using  predicted  current  inten¬ 
sity)  . 

LOCATION 

Sine  of  the  latitude  of  the  cyclone 
center  at  image  time. 

Sine  of  the  longitude  of  the  cyclone 
center  at  image  time. 

SEASON 

Julian  date  of  satellite  image  con¬ 
verted  to  a  sine  wave  and  phase 
shifted  to  day  80. 

DIURNAL  EFFECTS 

Image  time  converted  to  a  sine  wave 
and  phase  shifted  to  0600  LST  (local 
solar  time) . 

SPIRAL  GEOMETRY 

Inflow  angle  of  characteristic  spiral 
analysis  grid  of  the  cyclone. 
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characteristics  of  the  analysis  are  the  sum  of  the  variance 
explained  by  the  first  three  harmonics,  and  the  amplitude 
and  number  of  the  harmonic  explaining  the  maximum  amount  of 
variance.  As  the  explained  variance  spreads  over  the  spec¬ 
trum  the  spiral  bands  become  less  defined  over  the  analysis 
grid  and  the  signal  to  noise  ratio  is  reduced.  The  poten¬ 
tial  predictors  chosen  to  represent  the  spiral  signal  in 
the  analysis  are  the  number  of  harmonics  necessary  to  explain 
85%,  90%,  and  95%  of  the  observed  variance  of  the  IR  radiance 
field. 

Gentry  et  al.  (1980)  showed  that  measures  of  the  vari¬ 
ation  of  the  IR  radiance  field  were  important  for  predicting 
tropical  cyclone  intensity.  Therefore,  the  observed  vari¬ 
ance  and  the  standard  deviation  of  the  IR  radiance  field 
are  used  as  potential  predictors;  however,  they  are  calcu¬ 
lated  from  the  prewhitened  section  values.  To  include  a 
measure  of  the  variation  of  the  original  field  the  thres¬ 
hold  value  used  to  prewhiten  the  data  is  also  used  as  a 
potential  predictor. 

As  seen  in  Table  1,  the  contribution  of  the  results  of 
the  Fourier  analyses  to  the  field  of  potential  predictors 
has  been  further  reduced.  Of  the  potential  predictors  de¬ 
scribed  above,  it  was  arbitrarily  decided  to  use  only  their 
means  over  the  ten  sections  and  their  values  for  sections 
#2,  #5,  and  #8.  Although  this  decision  restricts  the 
analysis,  the  computations  are  done  more  efficiently  with 
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fewer  potential  predictors  and  the  interpretation  of  the 
results  is  less  complicated. 

3.2.2  Calculation  of  Regression  Equations 

The  stepwise  multiple  linear  regression  procedure  used 
to  relate  the  potential  predictors  to  tropical  cyclone  in¬ 
tensity  is  the  BMDP  program  P2R,  revision  6-81  (Dixon, 

1981) .  The  program  was  developed  at  the  Health  Sciences 
Computing  Facility,  UCLA,  under  NIH  Special  Research  Re¬ 
sources  Grant  RR-3.  For  all  equations  the  F-method,  forward 
stepping  algorithm  is  used  with  a  minimum  F-to-enter  value 
of  4.0,  minimum  F-to-remove  of  3.9,  and  a  tolerance  of 
0.05. 

In  the  present  study  no  variables  are  forced  into  the 
regression  equations;  all  equations  include  intercepts,  and 
no  limit  is  set  on  the  number  of  steps  taken.  The  number 
of  predictors  is  unrestricted.  Further  details  of  the 
calculation  of  multiple  linear  regression  equations  can  be 
found  in  Dixon  (1981)  and  Afifi  and  Azen  (1979)  . 

3.2.3  Independent  Tests 

The  linear  regression  equations  are  then  tested  on  an 
independent  data  set.  The  testing  method  used  is  to  com¬ 
pare  the  independent  data  set  to  its  estimate,  calculated 
using  the  regression  equations.  The  estimated  and  indepen¬ 
dent  data  sets  are  compared  using  the  BMDP  program  3D, 
revision  6-81  (Dixon,  1981)  .  BMDP  3D  calculates  the  pooled 
Student's  t-test  for  the  difference  between  the  means  of 
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two  groups,  the  Mahalanobis  distance  D2 ,  which  is  a  measure 
of  the  distance  between  two  groups,  and  Levene's  F-test  for 
equality  of  variances,  an  assumption  of  Student's  t-test. 
Probability  values  are  also  calculated. 

D  square  measures  the  distance  between  the  means  of 
two  groups.  It  is  defined  as 

D2  =  (Y  -  Y)TS-1(Y  -  Y)  (6) 

where  (Y  -  Y)T  is  the  transpose  of  the  matrix  difference  of 
the  means  of  the  independent  measurements  Y  and  the  results 

A 

of  the  regression  of  Y  on  X,  denoted  as  Y .  S  is  the  co- 
variance  matrix  of  the  two  groups  being  compared.  D  square 
can  be  shown  to  have  an  F  distribution  and  can,  therefore, 
be  assigned  a  probability  value  (p-value)  based  on  an  F 
distribution  table  with  one  degree  of  freedom  (Afifi  and 
Azen,  1979). 

Student's  t-test  measures  the  difference  between  the 
means  of  the  two  groups  and  can  also  be  assigned  a  p-value 
using  a  t  distribution  table  with  one  degree  of  freedom. 

A  critical  assumption  of  the  t-test  is  the  equality  of 
variances  of  the  two  groups.  This  assumption  is  tested  by 
Levene's  F-test  which  is  a  robust  test  performed  by  an 
analysis  of  variance  on  the  absolute  deviations  of  each  case 
from  its  cell  mean  (Brown  and  Forsythe,  1974) .  F  can  be 
shown  to  have  an  F  distribution  and  can,  therefore,  be 
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assigned  a  p-value  using  an  F  distribution  table  with  one 
degree  of  freedom  (Brown  and  Forsythe,  19  74)  . 
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4. 


DATA 


The  population  from  which  the  dependent  and  indepen¬ 
dent  data  sets  were  extracted  consists  of  97  digital  IR 
satellite  images  of  15  tropical  cyclones.  All  of  the  images 
are  2  n  mi  resolution  in  the  spectral  range  of  10.5-12.6  ym 
(thermal  IR  window) .  To  ensure  independent  data,  successive 
images  are  separated  by  at  least  12  h  (Gentry  et  al., 

1980) .  From  Table  2  we  see  that  six  tropical  cyclones  from 
the  1979  season,  three  East  Pacific  and  three  Atlantic,  for 
a  total  of  76  images,  comprise  the  dependent  data  set. 
Twenty-one  images  of  nine  tropical  cyclones  form  the  inde¬ 
pendent  data  set.  Table  3  shows  that  the  cyclones  range 
from  the  1979  to  the  1983  season.  Henri  (1979)  is  the  only 
Atlantic  cyclone  used  in  the  independent  data  set. 

In  Tables  2  and  3  the  image  times  are  divided  into 
0000  and  1200  GMT.  These  synoptic  times  are  only  approxi¬ 
mate  image  times,  used  to  show  the  independency  of  the  data. 

4.1  DEPENDENT  DATA  SET 

4.1.1  Satellite  Imagery 

The  digital  IR  satellite  data  for  76  images  processed 
with  TEDS  (Table  2)  were  acquired  from  GOES  data  tapes  pro¬ 
duced  on  the  Off-line  Data  Ingest  System  (ODIS)  by  the 
Space  Science  and  Engineering  Center  (SSEC) ,  University  of 
Wisconsin,  Madison.  These  ODIS  tapes  can  be  acquired 
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through  the  National  Environmental  Satellite  Data  and 
Information  Service  (NESDIS) ,  Washington,  D.C. 

The  Pacific  hurricanes  were  imaged  with  the  GOES-3 
satellite,  stationed  at  135  W  and  in  a  geosynchronous  equa¬ 
torial  orbit  at  an  altitude  of  approximately  35,800  km. 

GOES- 3  was  launched  16  June  1978  and  became  operational  as 
GOES-WEST  on  13  July  1978. 

The  IR  images  of  the  Atlantic  hurricanes  were  taken  by 
the  SMS-2  satellite  which  was  stationed  at  75  W  and  oper¬ 
ated  as  GOES-EAST.  SMS-2  was  also  a  geosynchronous  satel¬ 
lite,  capable  of  viewing  the  entire  disk  of  the  earth  with 
each  image.  The  radiometer  used  on  both  the  GOES-3  and 
SMS-2  satellites  was  the  Visible  Infrared  Spin  Scan  Radi¬ 
ometer  (VISSR)  which  provided  night/day  observations  of 
cloud  cover  and  Earth/cloud  radiance  temperature  measure¬ 
ments.  SMS-2  was  launched  6  February  1975  and  became  opera¬ 
tional  as  GOES-EAST  on  19  April  1979. 

4.1.2  Ground  Truth  Data 

Individual  descriptions  of  the  life  histories  of  the 
six  tropical  cyclones  comprising  the  dependent  data  set 
(Table  2)  are  given  in  Appendix  C.  The  positions  at  image 
times  of  Dolores,  Enrique,  and  Ignacio  are  supplied  by 
Gunther  (1980),  with  the  corresponding  positions  of  the 
Atlantic  hurricanes:  David,  Frederic,  and  Gloria  supplied 
by  Hebert  (1980) .  Because  some  of  the  images  used  in  this 
study  were  not  taken  at  0000  GMT  or  1200  GMT,  the  positions 
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are  interpolated  between  reports.  The  Akima  (1970)  inter¬ 
polation  method  is  used  because  it  keeps  the  first  deriva¬ 
tive  field  continuous  in  the  variable  being  analyzed,  and 
is  computationally  efficient. 

The  intensities  at  image  times  are  supplied  by  the 
Fleet  Numerical  Oceanography  Center's  archived  tropical 
cyclone  warning  data  file.  Instead  of  analyzing  the  inten¬ 
sity  of  the  non-standard  timed  images  individually  using 
the  Dvorak  technique,  the  Akima  interpolation  method  is 
again  used  and  rounded  to  the  whole  knot. 

4 . 2  INDEPENDENT  DATA  SET 
4.2.1  Satellite  Imagery 

The  tropical  cyclones  making  up  the  independent  data 
set  are  listed  in  Table  3.  The  Atlantic  hurricane  (Henri) 
was  imaged  by  the  SMS-2  satellite  and  Fefa  was  imaged  by 
GOES-3,  both  satellites  described  in  4.1.1.  The  images 
were  acquired  from  NESDIS .  The  IR  images  of  the  1982  season 
Pacific  hurricanes  were  also  taken  by  GOES-3.  The  IR 
images,  however,  were  collected  directly  on  NEPRF ' s  GOES 
acquisition  and  data  handling  system  (GADHS)  by  the  SPADS 
Department . 

The  digital  IR  satellite  images  for  the  two  1983  season 
Eastern  Pacific  tropical  cyclones  were  also  collected 
directly  by  the  SPADS  Department.  The  cyclones  were  imaged 
by  GOES-5,  stationed  at  75  W  and  operated  as  GOES-EAST. 
GOES-5  is  a  geosynchronous  satellite  but  uses  an  upgraded 
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radiometer:  VISSR  Atmospheric  Sounder  (VAS) .  However,  the 

upgrade  applies  to  the  sounding  capabilities  of  the  instru 
ment,  the  IR  images  are  still  sensed  with  the  VISSR.  GOES- 
5  was  launched  in  May  1981  and  became  operational  as  GOES- 
EAST  on  5  August  1981. 

4.2.2  Ground  Truth  Data 

The  position  at  image  times  of  Henri  and  Fefa  are 
interpolated  using  the  Akima  (1970)  method  from  Hebert 
(19  80)  and  Gunther  (19  80)  ,  respectively.  The  positions  and 
intensities  for  Adolf  and  Barbara,  and  the  intensities  for 
Henri  and  Fefa,  are  supplied  by  the  Fleet  Numerical  Ocean¬ 
ography  Center's  archived  tropical  cyclone  warning  data 
file.  The  intensities  at  image  times  are  interpolated  be¬ 
tween  reports  with  Akima 's  method.  The  positions  and  inten¬ 
sities  of  the  1982  hurricanes  are  supplied  by  Emil  Gunther 
of  the  Eastern  Pacific  Hurricane  Center,  Redwood  City,  and 
by  Hans  Rosendal,  Central  Hurricane  Center,  Honolulu,  both 
by  personal  communication.  The  data  provided  are  from 
post-analyzed  best- track  data  sets.  The  data  at  image 
times  are  again  interpolated  between  reports  using  Akima  s 
method. 

The  data  used  to  develop  the  12  and  24  h  intensity 
forecast  regression  equations  and  the  independent  data  used 
for  verification,  are  the  best-track  and  warning  intensi¬ 
ties  described  above.  Since  the  equations  are  neither 
developed  nor  verified  using  forecast  data  only  70  of  the 
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76  images  in  the  dependent  data  set  can  be  used  to  develop 
the  24  h  forecast  regression  equations.  The  last  image 
of  each  cyclone  has  to  be  eliminated  because  of  the  lack 
of  verification  data.  All  images  are  available  for  use  in 
the  0  and  12  h  forecast  equation  development.  Of  the  21 
independent  cases  Fefa,  24  August  1979,  showed  a  strikingly 
rapid  decay  (80  kt  over  36  h)  which  is  not  representative 
of  either  data  set  (mean  24  h  change  of  12  kt) .  For  this 
reason  Fefa,  24  August  1979,  has  been  eliminated  from  the 
results . 

Due  to  the  small  temporal  range  of  images  available 
for  each  tropical  cyclone,  individual  descriptions  are  not 
given.  Figure  11  is  a  plot  of  the  frequency  distribution 
for  the  current  intensity  (VMAX)  for  both  data  sets.  The 
histograms  are  divided  into  Dvorak  (1975)  classifications, 
the  T-number  and  associated  wind  speeds  being  the  abscissa. 
Note  the  ordinate  scale  change.  The  bimodal  distribution 
of  the  dependent  data  set  is  primarily  due  to  the  small 
sample  size;  however,  the  intensity  data  are  not  acquired 
from  direct  observations,  they  are  estimated  operationally 
using,  primarily,  the  Dvorak  (1975  and  1979)  technique. 
Therefore,  any  operational  bias  in  the  issuance  of  inten¬ 
sity  estimates  will  be  reflected  in  the  frequency  distribu¬ 
tion.  Figure  12  is  similar  to  Figure  11  except  that  is  is 
for  the  intensity  12  h  after  image  time  (VP12) .  The  plots 
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Figure  11.  Frequency  distributions  of  tropical  cyclone 
intensity  at  image  time  (VMAX)  for  the  dependent  (top) 
and  independent  (bottom)  data  sets.  Note  the  ordinate 
scale  change.  The  stratification  is  by  Dvorak  T-number. 
The  associated  wind  speed  is  in  knots. 
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Figure  12 .  Same  as  Fig.  11  except  for  intensity  12  h  after 
image  time  (VP12)  . 
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Figure  13 .  Same  as  Figs.  11  and  12  except  for  intensity 
24  h  after  image  time  (VP24) .  Only  70  cases  were  avail¬ 
able  in  the  dependent  data  set. 
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show  nearly  the  same  distribution  as  VMAX  except  for  an 
increase  in  the  number  of  storms  classified  T4  .  Since  most 
of  the  intensity  monitoring  techniques  of  2.2  rely  heavily 
on  persistence,  the  temporal  change  of  the  frequency  dis¬ 
tributions  are  slow,  as  evidenced  by  the  persistence  of  the 
bimodal  distribution  in  the  dependent  data  set.  Figure  13 
is  similar  to  Figures  11  and  12  except  for  the  change  of 
variable  to  intensity  24  h  after  image  time  (VP24) .  Figure 
13  shows  nearly  the  same  frequency  distribution  as  Figure 
12.  Note  that  only  70  cases  are  represented  by  the  depen¬ 
dent  data  set. 
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5.  RESULTS 


As  described  in  3.2.2,  regression  equations  were 
developed  to  specify  the  current  intensity  (VMAX) ,  and  to 
make  12  h  (VP12) ,  and  24  h  (VP24)  intensity  forecasts  for 
tropical  cyclones.  Following  Gentry  et  al.  (1980)  all  of 
the  cases  were  stratified  into  two  groups  depending  upon 
whether  or  not  the  cyclone  was  of  hurricane  intensity  12  h 
before  image  time  (VM12) :  VM12  £  65  kt  (LE65)  and 
VM12  >  65  kt  (GT65) .  This  stratification  also  ensured  the 
availability  of  at  least  23  cases  per  group  for  the  devel¬ 
opment  of  each  regression  equation.  The  regression  equa¬ 
tions  were  developed  for  the  three  dependent  variables 
using  both  the  stratified  and  the  unstratified  data;  there¬ 
fore,  the  model  consists  of  three  equations  for  each  depen¬ 
dent  variable.  One  equation  (ALL)  was  developed  using  all 
of  the  cases  and  can  always  be  used  to  specify  the  varia¬ 
ble.  One  of  the  two  remaining  equations  (LE65  and  GT65) 
can  be  used  to  again  specify  the  variable  depending  upon 
VM12  as  described  above. 

For  the  independent  testing  the  results  for  each 
dependent  variable  using  the  LE65  and  GT65  equations  were 
combined  to  form  the  group  GTLE65.  Therefore,  the  esti¬ 
mates  for  each  dependent  variable  form  two  groups;  the 
group  ALL,  calculated  using  the  ALL  equation,  and  the 
group  GTLE65,  calculated  using  the  LE65  and  GT65  equations. 
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The  independent  test  cases  were  also  divided  into  the 
groups  ALL,  LE65,  and  GT65,  and  the  recombination  was  per¬ 
formed  the  same  as  for  the  dependent  data  set.  The  strati¬ 
fication  for  the  20  case  independent  data  set  is  nine  cases 
LE65  and  eleven  cases  GT65. 

Table  4  is  a  list  of  the  descriptions  of  the  predictors 
picked  by  BMDP2R  for  the  linear  regression  equations  and 
the  descriptions  of  the  dependent  variables.  The  predictors 
were  selected  from  the  field  of  potential  predictors  dis¬ 
cussed  in  3.2.1.  Tables  5a  and  5b  are  descriptive  statis¬ 
tics  for  the  predictors  of  Table  4.  Table  5a  represents 
the  dependent  data  set  and  Table  5b  represents  the  indepen¬ 
dent  data  set. 

Comparing  the  two  data  sets,  we  find  from  the  differ¬ 
ences  in  the  kurtosis  that  the  dependent  data  has  a  higher 
probability  of  extreme  values  than  the  independent  data 
set.  The  means  of  the  variables  agree  quite  closely  ex¬ 
cept  for  the  measurements  of  persistence  (DVM12  and  DV24) 
and  ALFA.  This  implies  the  independent  data  set  is  com¬ 
prised  of  tropical  cyclones  that,  on  the  average,  devel¬ 
oped  more  rapidly  than  those  making  up  the  dependent  data 
set.  The  mean  of  ALFA  is  7.5°  less  for  the  more  rapidly 
developing  cyclones. 

Comparing  ALFA  to  Swadley's  (1983)  results,  using  the 
polar-log  spiral  equation  but  tracing  spiral  cloud  bands, 
we  find  that  the  mean  inflow  angles  (~28°)  are  nearly 
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Table  4.  Description  of  predictors  selected  by  BMDP2R  for 
the  multiple  linear  regression  equations  and  descriptions 
of  the  dependent  variables. 


NAME_ DESCRIPTION 

X90M  number  of  harmonics  necessary  to  explain  90%  of 

the  observed  variance  for  the  mean  of  all  sections 

V32  percentage  of  the  variance  explained  by  the  first 

three  harmonics  of  section  2 

HEMVM  number  of  the  harmonic  explaining  the  maximum 
variance  for  the  mean  of  all  sections 

AMP 8  amplitude  (IR  counts)  of  the  harmonic  explaining 

the  maximum  variance  for  section  8 

AMPM  same  as  AMP8  except  for  the  mean  of  all  sections 

2 

TVAR5  total  observed  variance  (IR  counts)  for  section  5 

TVAR8  same  as  TVAR5  except  for  section  8 

SD2  standard  deviation  (IR  counts)  of  section  2 

THLD8  threshold  value  (IR  counts)  of  section  8 

TIMES  image  time  normalized  to  0600  LST  and  converted 
to  a  sine  wave;  TIMES  =  sin  (  (  2tt/24)  T)  where 
T  =  MOD (LT-6+24 , 24)  and  LT  is  the  local  solar 
time  of  the  image 

VM12  intensity  (kt)  12  h  before  image  time 

DVM12  change  in  intensity  (kt)  between  24  and  12  h  before 

image  time;  DVM12  =  VM12  -  VM24  where  VM24  is  the 
intensity  24  h  before  image  time 

DV24  change  in  intensity  (kt)  24  h  before  image  time 

using  VMAX  estimated  by  the  appropriate  regression 
equation;  DV24  =  VMAX  -  VM24 

ALFA  inflow  angle  (degrees)  of  the  characteristic  spiral 


VMAX 
VP  2  4 
VP12 


intensity 

(kt) 

at 

intensity 

(kt) 

24 

intensity 

(kt) 

12 

image  time 
h  after  image  time 
h  after  image  time 


Dependent 

Variables 
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Table  5a.  The  mean,  standard  deviation,  skewness,  and 
kurtosis  for  all  of  the  predictors  picked  by  BMDP2R 
and  the  dependent  variables.  The  variables  are  des¬ 
cribed  in  Table  4. 


VARIABLE 

NAME 

MEAN 

STANDARD 

DEVI  AT  ION 

SKEWNESS 

KU&T0SIS 

X90M 

•  6.447 

1 . 670 

.791 

.219 

V  32 

75.907 

11.635 

-.858 

-.188 

V3M 

79.444 

5. 8B7 

-.752 

.314 

HEMVM 

1.195 

.  239 

1.028 

-.164 

AMP8 

211.964 

16.039 

-2.467 

8.436 

A  MPM 

309.362 

15.465 

-1.441 

2.080 

TVAR5 

289.140 

374.623 

3.390 

13.335 

T  VAR  8 

255.034 

340.741 

2.444 

6.357 

S  02 

24.525 

11. 487 

.267 

-.737 

THLD8 

171.289 

30.092 

-.737 

-.226 

TIMFS 

.019 

.  42  1 

-.146 

-1.361 

VM1  2 

65.395 

29.344 

.641 

-.273 

DVM12 

4.908 

11.466 

.472 

1.403 

DV24 

9.210 

11.251 

-1.454 

6  •  445 

ALFA 

28. 364 

9.  B88 

.296 

.085 

V  M  A  X 

69.487 

28.183 

.634 

-.217 

VP  2  4 

75.286 

27.467 

.708 

-.401 

V  P 1  2 

72.329 

27.656 

.681 

-.331 

Table 

data 

5b.  Same  as 
set. 

Table  5a  except 

for  the  independent 

VARIA3LE 

NAME 

mean 

S  TANDARD 

DE  VIATIQN 

SKEWNESS 

KURTOSIS 

X90M 

6.935 

1.661 

1.213 

1.658 

V  3  2 

7  2 . 67  6 

1 3.266 

-.534 

.028 

V3M 

76.494 

7.6*>9 

-1.156 

1.351 

HEMVM 

i.  n 

.3  24 

1.403 

.422 

AMPP 

2  j  3.46  fir 

8.210 

.448 

-.039 

A  M  p  M 

210. 7 0C 

5.9  or 

- .  664 

-.956 

TVAR5 

259.179 

264.633 

1.294 

.938 

T  V  A  P  8 

i  n.  600 

19  4,435 

2.003 

3.587 

S  0  2 

20.236 

7.164 

.400 

-.773 

THLD8 

177.900 

19.378 

-.491 

-.  528 

TIMES 

.  59  3 

.453 

-1.417 

.845 

VM12 

67.250 

23.420 

.224 

-.492 

DVM12 

8.3  00 

12.215 

.136 

-.958 

DV24 

20.114 

9.325 

-.673 

-1.045 

ALFA 

21.926 

4.327 

.300 

-.656 

V  M  A  X 

73.650 

2  1.176 

.321 

-.906 

VP24 

68.250 

18.487 

-.222 

-.813 

V  P 1 2 

72.800 

1 5.917 

-.202 

-.728 
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equal.  Therefore,  the  polar-log  spiral  representation  of 
tropical  cyclone  cloud  bands  appears  to  be  valid.  This  is 
expected  since  local  polar  coordinates  are  a  good  approxi¬ 
mation  to  local  spherical  coordinates  over  the  distance 
scales  of  a  tropical  cyclone. 

5.1  CURRENT  INTENSITY  (VMAX) 

Results  for  current  intensity  (VMAX)  from  the  stepwise 
multiple  linear  regression  program  BMDP2R  and  the  results 
from  the  two  groups  comparison  program  BMDP3D  are  shown  in 
Table  6a  (ALL)  and  Table  6b  (GTLE65) .  In  these,  and  simi¬ 
lar,  tables  the  percentage  of  the  variance  explained  by 
the  regression  equation  (PVE)  is  the  multiple  correlation 
coefficient,  R  ,  squared  and  multiplied  by  one-hundred. 

The  F  ratio  included  in  the  analysis  of  variance  (ANOVA)  is 
a  test  of  the  significance  of  the  coefficients  of  the  pre¬ 
dictors  in  the  regression  equations.  The  potential  predic¬ 
tors  are  selected  in  such  a  manner  as  to  maximize  it.  In 
general,  the  coefficients  in  the  regression  equations  are 
significant  at  the  5%  level  and  are  listed  in  the  equations 
from  left  to  right  in  order  of  their  relative  significance. 

The  equations  for  ALL  and  GT65  show  that  persistence  is 
a  more  important  predictor  for  current  intensity  than  are 
the  spiral  parameters;  the  average  contribution  to  VMAX 
(ALL)  from  the  predictor  SD2  is  a  decrease  in  intensity  of 
5.5  kt  and  the  GT65  equation  contains  only  persistence. 

The  estimation  of  the  current  intensity  for  tropical 
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Table  6a.  Results  of  the  regression,  ANOVA,  and  comparison 
of  two  groups  analyses  for  the  group  VMAX  ALL. 


EQUATION 


VMa”~(ALL)  «  21.7817?  +  0.7911*VM12  +  0 . 5 45 1 7*0 VM 1 2  -  C.27336+S02 


MULTIPLE  P 
MULTIPLE  R-SQUARE 
STD.  ERROR  OF  EST. 


.97966 

.95974 

6.66785 


ANALYSIS  OF  VARIANCE 


SUM 

REGRESSION 

RESIDUAL 

OF  SQUARES 
3968.550 
711.364 

DF 

3.0 

16.0 

mean  square 

2989.517 

44.460 

F  RATIO 
67.240 

STATISTICS 

MAHALANOBIS  D  SQUARE 

T  (POOLED) 

F  (FOR  VARIANCES)  LEVENE 

P-VALUE 
.000?  .9618 

-.05  .9618 

.13  .7229 

DF 

It  38 
It  38 
It  38 

Table  6b.  Same  as  Table  6a  except  for  the  group  VMAX 
GTLE65. 


EQUATION 

VMAx"7lE66)  =  -?7. 12597  +  0.67355*VM1’  ♦  0  .  ?? 71  A « *  6. 7213 6* TIMES  -  0 . 3409  ?  *  AL  F  A 
VM4X  (GT65)  *  31.63646  *  C  .  4  ?  706*  V  *  1  2  +  0 . 7 1 402* 0 V 'U  2 


MULT  I  PL  F  P 
MULTIPLE  P-SO'JARF 
STD.  FRR09  OF  EST 


.94339 

.89092 

11.71577 


ANALYSIS  OF  VARTANCF 


REGRESSION 

RESIDUAL 

STATISTICS 


SUM  OF  SQUARES 
8968.550 
2196.150 


OF 


MAHALANOPIS  0  SOUARF 
T  (POOLED) 

F  (FOR  VARIANCES)  LFVENE 


^  P  A  N  S  OU  A  p  F 

F  RATIO 

0 

?9«9.517 

21.780 

0 

137.259 

P-VALUF 

DF 

0096 

.  7590 

It  ?8 

31 

.7590 

It  38 

03 

.  8707 

It  38 
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cyclones  below  hurricane  strength  (LE65)  depends  on  the 
cyclone's  spiral  geometry  derived  from  the  cloud  bands, 
persistence,  and  image  time.  The  inclusion  of  TIMES  in 
the  VMAX  (LE65)  equation  suggests  a  diurnal  variation  in 
the  current  intensity  of  the  weak  storms.  Browner  et  al . 
(1977)  investigated  the  diurnal  oscillation  of  the  cirrus 
canopies  of  four  tropical  cyclones  using  satellite  derived 
IR  data,  and  found  a  negative  correlation  between  storm 
intensity  and  the  amplitude  of  the  canopy  oscillations. 

Thus,  it  appears  that  for  weak  storms  the  intensity  and 
cloud  top  temperature  are  related  on  a  diurnal  scale. 

The  equations  also  show  that  a  combination  of  persist¬ 
ence  and  spiral  parameters  is  capable  of  explaining  a  large 
percentage  of  the  variance  of  VMAX.  The  maximum  PVE  is  95% 
for  ALL  and  the  standard  error  (SDE)  is  6  kt.  The  PVE  for 
GTLE 6 5  is  89%  and  the  SDE  is  12  kt.  The  PVE  for  both  ALL 
and  GTLE65  is  larger  than  the  76%  reported  by  Swadley 
(1983)  for  a  19  case  dependent  data  set  and  using  no  per¬ 
sistence.  Swadley  presented  no  independent  results.  Both 
equations  exceed  the  noise  level  of  Dvorak  current  inten¬ 
sity  estimates  as  reported  by  Erickson  (1972)  and  Sheets 
and  Grieman  (1975) .  Erickson  verified  Dvorak  estimates 
using  aircraft  reconnaissance  data  for  300  cases  and  re¬ 
ported  a  PVE  of  78%  and  an  RMS  error  (RMSE)  of  15  kt. 

Sheets  and  Grieman  verified  480  Dvorak  estimates  against 
the  corresponding  best-track  data  and  found  that  92%  of  the 


52 


estimates  were  correct  to  within  one  T-number  and  that  the 


internal  consistency  of  the  technique  between  experienced 
analysts  was  1.5  T-numbers. 

Table  7  is  a  tabulation  of  all  the  independent  re¬ 
sults.  Along  the  bottom  are  the  means,  standard  deviations 
(STDV) ,  total  variances  (TVAR) ,  RMS  errors,  and  PVE  for  the 
groups.  The  RMSE  for  ALL  and  for  GTLE65  are  both  less  than 
the  RMSE  reported  by  Erickson  and  the  RSME  for  ALL  is  four 
knots  smaller  than  Swadley's  (1983)  results.  Figure  14  is 
a  plot  of  the  residuals  of  Table  7  against  the  estimated 
VMAX .  The  grouping  in  the  higher  wind  speeds  is  an  arti¬ 
fact  of  the  Dvorak  stratification  in  persistence  for  the 
predictors  in  the  GT65  equation. 

Returning  to  Tables  6a  and  6b  we  find  that  the  D  square 
distances  between  the  means  of  the  estimates  and  observa¬ 
tions  are  very  small,  and  that  the  t-test  for  GTLE65  is  six 
times  that  for  ALL.  The  p- values  confirm  that  we  must  ac¬ 
cept  the  null  hypothesis  (Hq)  for  both  groups.  In  other 
words,  the  probability  is  96%  for  ALL  and  76%  for  GTLE65 
that  the  two  groups  are  members  of  the  same  population  as 
VMAX.  Although  the  Levene  F-test  has  smaller  p-values,  we 

still  cannot  reject  H  . 

o 

From  an  F  distribution  table,  we  find  that  the  F 
ratios  for  both  ALL  and  GTLE65  are  significant  at  the  1% 
level.  This  means  that  there  is  only  one  chance  in 


53 


Table  7.  Tabulated  results  of  the  regression  analyses  for  VMAX,  VP12,  and 
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Figure  14 .  Residuals  plotted  against  estimated  VMAX  ALL 
(top  graph)  and  estimated  VMAX  GTLE65  (bottom  graph) . 
The  residuals  are  given  by  Y-Y. 
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one  hundred  that  a  random  set  of  predictors  could  produce 
these  results. 

In  summary,  the  results  for  VMAX  show  that  a  combina¬ 
tion  of  persistence  and  spiral  parameters  is  statistically 
significant  in  correctly  estimating  the  current  intensity 
of  tropical  cyclones.  The  errors  in  the  model  are  at  the 
noise  level  of  the  verification  data,  and  are  better  than 
other  published  techniques. 

5.2  12  h  FORECAST  INTENSITY  (VP12) 

Tables  8a  and  8b  show  the  results  of  BMDP2R  and  BMDP3D 
for  VP12.  Note  VMAX  in  the  regression  equations  are  those 
estimated  by  the  corresponding  equations  and  groups  for 
consistency.  The  results  for  the  12  h  forecast  data  are 
not  as  impressive  as  those  for  the  current  intensity.  The 
PVE  for  ALL  (60%)  is  twice  that  for  GTLE65.  The  standard 
errors  are  at  the  class  limits  for  the  larger  T-numbers  of 
the  Dvorak  classification  scheme.  This  result  is  also 
indicated  graphically  from  the  contingency  tables  of  Table 
9.  Considering  total  errors,  GTLE65  performs  better  than 
ALL,  ten  errors  to  twelve,  however,  the  number  of  two  class 
errors  is  larger  for  GTLE65  (4)  than  for  ALL  (1) ,  which 
means  VP12  was  correctly  forecast  by  ALL  19  out  of  20  times 
within  one  T-number  category. 

A  skill  score  S  is  defined  as 
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Table  8a.  Same  as  Table  6a  except  for  the  group  VP12  ALL 


EQUATION 

Vpi2  (ALL)  *  -152.49974  +  0.75016+VMAX  +  0 . 67 8 9 3* AN DM  +  2 1 . 1 4 516* HE M VM  ♦  0. 0 1 8  36* T  V  AR  5 


MULTIPLE  R 

.77908 

MULTIPLE  R-SQUARE 

.60697 

STD.  ERROR  OF  EST. 

16.42223 

ANALYSIS  OF  VARIANCE 

SUM 

OF  SOUARFS  OP 

MEAN  SQUARE 

F  RATIO 

REGRESSION 

5067.200  4. 

0 

1266. 800 

4.697 

RESIDUAL 

4045.346  15. 

0 

269.690 

STATISTICS 

o-VALUE 

OF 

MAHALANOBIS  D  SQUARE 

• 

0691 

.4110 

If  38 

T  (POOLED) 

—  • 

83 

.4110 

If  38 

F  (FOR  VARIANCES)  LFVENE 

• 

13 

.7174 

If  38 

Table  8b.  Same  as  Table  6a  except  for  the  group  VP12 
GTLE65. 

EQUATION 

VP12  (LE65)  =  -130. 3^8)9  +  0 . 99 5 2 5* VIA  X  +  0.52087*AMP“  +  2 1 . 2 6 3 24*HE M VM  +  0 • 0 1 3 5 5*T V AR 5 
VP12  (GT65)  *  -217.34952  +  0 . 5  P'2  26  *  VM  A  X  +  C..07099*TVAR8  +  1 , 36664+THLD8 


MULTIPLE  P 
MULTIPLE  R-SQUARF 
STO.  ERROR  OF  FST . 


.55315 
.  30597 
23.09397 


ANALYSIS  OF  VARIANCE 


REGRESSION 

RESIDUAL 

STATISTICS 


MAHALANOBIS  D  SQUARE 
T  (POOLED) 

F  (FOR  VARIANCES)  LEVE 


SUM  OF  SOUARFS  OF  MEAN  SQUARE 

F  RATIO 

5067.200  3.5 

1447.771 

2.715 

8266.639  15.5 

533.33? 

P-VALUE 

OF 

.0650 

.  4252 

1,  38 

-.81 

.425? 

1, 38 

NP  .45 

.  5073 

1>  38 
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Table  9.  Contingency  tables  for  VP12  ALL  (left)  and  VP12  GTLE65  (right)  classified 
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where  P  is  the  number  of  correct  forecasts  by  persistence, 

C  is  the  number  of  correct  forecasts  by  the  model,  and  T  is 
the  total  number  of  forecasts.  The  S  for  ALL  is  -0.33  and 
the  S  for  GTLE65  is  -0.11.  These  results  indicate  the 
models  are  not  skillful  in  predicting  VP12  within  the  cor¬ 
rect  Dvorak  T-number  category  when  compared  to  persistence. 
These  results  cannot  be  compared  because  there  are  no  pub¬ 
lished  12  h  forecast  tropical  cyclone  intensity  error 
statistics . 

From  Table  7  we  find  that  the  RMSE  for  ALL  is  14  kt, 
one-third  better  than  the  RMSE  for  GTLE65 .  This  is  not  re¬ 
flected  in  the  skill  score.  The  residuals  are  plotted 
against  the  predicted  VP12  in  Figure  15.  The  results  show 
a  randomness  that  is  usually  associated  with  an  adequate 
model  (Afifi  and  Azen,  1979) . 

From  Table  8a  we  find  that  the  F  ratio  for  ALL  is  sig¬ 
nificant  at  the  5%  level.  The  F  ratios  are  an  order  of 
magnitude  less  than  those  for  VMAX,  but  the  p-values  still 
require  that  we  accept  Hq;  the  two  groups  are  most  likely 
from  the  same  population  as  VP12. 

From  Table  7  we  can  see  that  the  difference  between 
the  means  is  larger  for  GTLE65  than  it  is  for  ALL.  This  4 
kt  bias  in  the  model  is  due  to  the  difference  between  the 
mean  rate  of  intensity  change  between  the  dependent  and 
independent  data  sets,  discussed  earlier  in  this  chapter, 


59 


£ 


CO 

u 


CO 

£ 


Figure  15.  Same  as  Fig.  14  except  for  the  population  VP12 
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and  should  only  be  removed  after  testing  the  model  on  a 
larger  independent  data  set . 

It  is  interesting  to  note  that  the  form  of  the  equa¬ 
tions  for  ALL  and  LE65  is  the  same;  this  is  carried  through 
to  VP24  (ALL),  discussed  next.  This  reflects  the  consis¬ 
tency  of  the  spiral  representation  of  the  banding  structure 
of  the  tropical  cyclone  and  validates  the  credibility  of 
the  model . 

5.3  24  h  FORECAST  INTENSITY  (VP24) 

The  results  of  the  regression  analyses  and  the  indepen¬ 
dent  test  for  VP24  are  shown  in  Tables  10a  and  10b.  The 
form  of  the  regression  equation  for  ALL  was  mentioned  pre¬ 
viously.  A  reason  that  this  form  of  equation  persists  is 
that  the  Fourier  analyses  always  include  the  coldest  cloud 
tops.  Gentry  et  al.  (1980)  showed  that  a  333  km  radial  ring 
about  the  storm  center  enclosed  the  area  of  coldest  cloud 
tops,  the  overshooting  convection  in  and  around  the  eye 
wall.  Gentry  et  al .  also  showed  the  average  cloud  top 
temperature  for  the  radial  ring  to  be  a  significant  predic¬ 
tor  for  24  h  intensity  forecasts.  This  importance  is  indi¬ 
cated  in  the  present  study  by  the  persistence  in  the 
equations  of  the  number  and  amplitude  of  the  harmonic  ex¬ 
plaining  the  maximum  variance  for  the  mean  of  all  sections. 
These  predictors  suggest  the  importance  of  the  high  IR 
values  of  the  coldest  cloud  tops  as  intensity  predictors. 
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Table  10a.  Same  as  Table  6a  except  for  VP24  ALL 


EQUATION 

VP24  (ALL)  *  -204.98015  ♦  0.48934+VMAX  +  0.96023+AMPM  ♦  3 2 . 69 795*HE M VM  +  0 . 01992+TV AR 5 


MULTIPLE  R 
MULTIPLE  P-SQUARE 
STD.  ERROR  OF  EST. 


.44945 
.20201 
26. 59049 


ANALYSIS  OF  VARIANCE 


REGRESSION 

RESIDUAL 

STATISTICS 


SUM  OF  SQUARES 
6835.750 
10605.313 


MAHALAN08IS  D  SQUARE 
T  (POOLED) 

F  (FOR  VARIANCES)  LEVENE 


OF  MF  AN  SQUARE 

F  RATIO 

4.0 

1708.938 

2.417 

15.0 

707.054 

P-VALUE 

OF 

.4195 

.0475 

If  38 

-2.05 

.0475 

1,38 

2.  59 

.1159 

1,38 

Table  10b.  Same  as  Table  6a  except  for  VP24  GTLE65 . 


EQUATION 

VE?4  (LE65)  *  -10. 71’ 31  +  0.72176  +  VMAX  +  0.7412*0V24+  0.47°04*V'52 
VP2A  ( G T 6 5  )  *  23.51729  +  2 . 5 46 1 8  *  A -  1 0 . 3 a 2 06+ X 90M  -  5.17114*V3U 


WULTTPLF  c  .51542 
MULTIPLF  P-SQUARE  .26566 
STD.  ERROR  OF  EST.  37.63501 


ANALYSIS  nF  VARIANCE 


REGRESSION 

RESIDUAL 

SUM  OF  SQUARES 
6335.750 
22662.299 

OF  m 

3.0 

16.0 

PAN  SQUAPF 
2278.583 
1416.394 

F  RATIO 
1. 609 

STATISTIC* 

MAHALANOPIS  D  SQUARE 

T  (POOLED) 

F  (FOR  VARIANCES)  LEVENE 

.  5541 
-2.35 

0.  00 

p-VALUE 

.0238 

.0238 

.9995 

DF 

Ip  38 
1,  38 
1#  38 
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This  agrees  with  a  conclusion  of  Swadley  (1983),  who  found 
the  maximum  IR  value  of  the  coldest  cloud  top  to  be  a  sig¬ 
nificant  predictor  of  tropical  cyclone  current  intensity. 

The  equation  for  GT65  is  interesting  because  it  is 
composed  only  of  spiral  parameters.  Persistence  is  found 
not  to  be  a  significant  predictor  for  24  h  forecasts  of 
intense  tropical  cyclone  maximum  wind  speed.  This  agrees 
with  the  results  of  Sheets  and  Grieman  (1975)  who  found  that 
the  Dvorak  forecasting  method  (which  relies  highly  on  per¬ 
sistence)  works  worst  for  intense  storms. 

The  PVE  for  ALL  is  20%  and  for  GTLE65  is  26%,  and  is 
much  less  than  that  explained  at  12  h.  From  Table  11  we 
find  that  the  number  of  class  errors  larger  than  one  is 
greater  for  GTLE65,  indicating  a  better  forecast  by  ALL 
when  compared  to  persistence,  even  though  ALL  has  a  lower 
PVE.  This  result  is  also  evident  in  the  skill  score,  which 
is  0.08  for  ALL  and  -0.33  for  GTLE65.  The  prediction  for 
VP24  made  by  using  the  ALL  equation  is  the  only  forecast 
more  skillful  than  persistence. 

From  Tables  10a  and  10b,  we  see  that  the  SDE  for  ALL  is 
26  kt  and  the  SDE  for  GTLE65  is  37  kt,  much  greater  than  one 
Dvorak  T-number  interval.  The  statistics  for  the  analysis 
of  the  two  groups  show  that  the  p-values  for  the  pooled 
Student's  t-test  and  D  square  are  low  enough  to  reject  Hq. 

At  the  4%  level  for  ALL  and  at  the  2%  level  for  GTLE65  we 
must  peremptorily  assume  that  the  two  groups  are  not  from 
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the  same  population  as  VP24 .  The  results  indicate  that 
these  regression  equations  have  no  usefulness  as  forecast 
aids . 

In  a  verification  of  the  Dvorak  (1975)  forecast  method 
against  aircraft  reconnaissance  data,  Erickson  (1972)  found 
that  the  24  h  forecast  technique  explained  31%  of  the  vari¬ 
ance  over  178  observations.  This  implies  that  the  noise 
level  of  the  Dvorak  forecasting  method  is  near  the  skill 
level  of  the  objective  forecast  method  presented.  These 
results  are  similar  to  those  of  Sheets  and  Grieman  (1975) 
who  found  that  the  24  h  forecast  errors  of  the  Dvorak  method 
were  of  the  same  magnitude  as  the  errors  made  by  forecasts 
using  climatology.  Shewchuk  and  Weir  (1980)  showed  that  the 
mean  error  for  the  Dvorak  (1979)  forecasting  method  over 
400  cases  was  8  kt  when  verified  against  best-track  data. 
Their  results  are  similar  to  the  errors  for  24  h  intensity 
forecasts  made  by  using  the  objective  technique  described 
by  Gentry  et  al.  (1980) .  Their  24  h  forecast  mean  absolute 
error  was  reported  as  14  kt,  which,  from  Table  7,  compares 
to  RMS  errors  of  23  kt  for  ALL  and  33  kt  for  GTLE65 .  The 
results  presented  above  indicate  that  the  24  h  objective 
forecast  technique  presented  in  this  study  is  not  as  skill¬ 
ful  as  the  other  satellite-based  tropical  cyclone  forecast¬ 
ing  techniques  in  the  literature. 

We  can  also  see  from  Table  7  that  ALL  can  only  repro¬ 
duce  about  one-half  of  the  total  variance  observed,  while 
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Figure  16 . 


Same  as  Fig. 


14  except  for  the  population  VP24. 
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the  TVAR  of  GTLE65  is  very  close  to  that  observed  for  the 
tropical  cyclones.  The  poor  results  shown  by  PVE  and  RMSE 
are  partly  due  to  the  differences  between  the  means  of  the 
groups.  The  bias  in  ALL  of  +10  kt  and  in  GTLE65  of  +14  kt 
are  close  to  the  bias  of  +11  kt  Gentry  et  al.  (1980)  found 
in  their  results.  The  biases  are  evident  in  the  large  num¬ 
ber  of  negative  residuals,  graphed  in  Figure  16  against  the 
forecasts  for  VP24.  The  bias  is  particularly  evident  for 
the  large  values  in  GTLE65 .  The  bias  in  the  forecast  is 
too  large  to  be  attributed  solely  to  differences  between 
the  data  sets  as  it  was  for  VP12.  The  results  indicate 
that  other  factors  relating  to  tropical  cyclone  intensity 
have  to  be  taken  into  account  in  order  to  develop  better 
24  h  predictive  techniques  (Afifi  and  Azen,  1979) .  With  a 
limited  data  set  further  conclusions  about  the  cause  of  the 
bias  are  not  possible. 
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6 .  CONCLUSIONS  AND  RECOMMENDATIONS 


6.1  CONCLUSIONS 

The  objective  of  this  study  was  to  determine  a  quan¬ 
titative  relationship  between  satellite  measurements  of 
the  IR  radiance  field  and  the  intensity  of  tropical  cy¬ 
clones.  A  spiral  analysis  technique  was  used  to  quantify 
the  spiral  patterns  of  the  cloud  bands  in  the  IR  data  and 
multiple  linear  regression  techniques  were  used  to  relate 
the  derived  spiral  parameters  to  tropical  cyclone  inten¬ 
sity.  The  results  presented  validate  the  hypothesis  that 
the  spiral  structure  of  cloud  bands  associated  with  tropi¬ 
cal  cyclones  can  be  characterized  by  fitting  spherical-log 
spirals  to  the  IR  radiance  field  of  their  cloud  tops;  how- 
ever ,  the  polar— log  spiral  representation  of  cloud  spirals 
used  by  other  researchers  appears  to  be  valid  for  tropical 
cyclones  located  near  the  satellite  sub— point. 

The  regression  equations  developed  show  significant 
skill  in  estimating  current  tropical  cyclone  intensity  and 
in  making  12  h  intensity  forecasts.  This  nowcasting  skill 
is  unique  because  the  technique  presented  is  a  stand-alone 
objective  analysis  and  forecast  method.  The  only  input  is 
a  current  IR  image  and  the  tropical  cyclone's  intensity 
for  the  previous  12  and  24  h. 

There  is  a  positive  bias  in  the  forecasted  intensities 
that  is  similar  to  the  bias  found  by  Gentry  et  al.  (1980) . 
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This  suggests  that  other  variables  not  considered  here  are 
necessary  to  correctly  predict  tropical  cyclone  intensity 
or  that  the  IR  radiance  field's  signal  to  noise  ratio  is 
too  low  for  the  signal  to  be  properly  detected. 

The  predictors  that  occur  most  frequently  in  the  re¬ 
gression  equations  are  related  to  the  maximum  IR  count 
found  in  the  central  dense  overcast  area  around  the  center 
of  tropical  cyclones.  This  result  was  also  found  by  Swadley 
(1983)  to  be  important  for  current  intensity  estimates,  and 
by  Gentry  et  al.  (1980)  to  be  important  for  intensity 
forecasts . 

A  negative  result  regarding  predictive  skill  was  found 
for  the  variable  ALFA,  the  inflow  angle  of  the  character¬ 
istic  spherical-log  spiral  analysis  grid  for  the  tropical 
cyclone.  This  supports  the  findings  of  Peterson  et  al . 
(1983)  and  contradicts  the  findings  of  Lahiri  (1981)  and 
Swadley  (1983) . 

6.2  RECOMMENDATIONS  FOR  FUTURE  WORK 

The  results  presented  above  demonstrate  that  more  re¬ 
search  needs  to  be  undertaken  in  the  area  of  satellite  de¬ 
rived  tropical  cyclone  intensity  estimates  and  forecasts. 

A  natural  expansion  of  this  study  would  include  other 
large  scale  potential  predictors  available  on  an  IR  image. 
Some  examples  of  additional  potential  predictors  are  the 
size  of  the  cyclone  measured  along  and  perpendicular  to  a 
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set  of  predetermined  axes;  the  areal  extent  of  a  certain 
equivalent  blackbody  cloud  top  temperature  corrected  for 
diurnal  oscillations  using  the  results  of  Browner ' s  et  al. 
(1977)  study;  the  width  and  azimuthal  extent  of  large  cloud 
bands  using  methods  developed  by  Dvorak  (1975,  1979);  and 
sea  surface  temperatures  in  the  cloud  free  areas  around 
the  cyclone  which  are  important  for  longer  range  forecasts. 
In  addition,  one  could  include  more  of  the  cloud  covered 
area  in  the  analysis  grid.  A  typical  spiral  grid  used  in 
this  study  covered  only  about  one-third  of  the  cloud  field 
of  the  tropical  cyclone. 

A  fitted  spiral  with  a  variable  offset  radius  could  be 
added  to  TEDS  in  order  to  more  accurately  determine  the 
spiral  characteristics  of  the  cloud  bands.  Energy  consider¬ 
ations  constrain  the  radial  inflow  so  the  inflowing  air 
cannot  penetrate  to  the  geometric  center  of  the  storm.  At 
a  typical  distance  of  10  to  100  km,  the  air  turns  upward  in 
a  ring  of  convection  called  the  eye  wall;  and  here  is  where 
the  strongest  wind  speeds  typically  occur  (Anthes,  19  82)  . 
This  radius  of  maximum  wind  speed  is  an  important  parameter 
to  be  able  to  predict,  and  its  correlation  with  the  offset 
radius  should  be  investigated.  In  addition,  the  behavior 
of  the  spiral  energy  spectrum  should  be  investigated. 

This  study  concentrates  on  information  available  from 
IR  images  of  tropical  cyclones,  consequently,  the  cloud 
bands  analyzed  are  composed  of  the  coldest,  upper-level 
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clouds  discernible  on  an  image.  But  the  widely  used  Dvorak 
(1975)  method  for  estimating  tropical  cyclone  intensity 
relies  on  visible  imagery  for  the  spiral  analysis  of  low- 
level  cloud  bands.  TEDS  also  has  the  capability  of  fitting 
spherical-log  spirals  to  cloud  bands  on  visible  images;  IR 
images  are  only  necessary  for  the  quantitative  data  to  be 
Fourier  decomposed.  This  capability  could  be  exploited  by 
analyzing  image  pairs:  visible  and  IR  images  temporally 
spaced  no  further  than  thirty  minutes  apart.  This  idea 
could  be  expanded  further  by  analyzing  image  sets  and  in¬ 
cluding  the  temporal  changes  of  the  spiral  parameters  as 
potential  predictors. 

Finally,  both  the  dependent  and  independent  data  sets 
should  be  expanded  from  their  sizes  in  the  present  study. 
More  significance  could,  therefore,  be  placed  on  the 
statistical  arguments. 
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APPENDIX  A 


TEDS  SOFTWARE  DOCUMENTATION 

Figure  Al  is  the  functional  flow  diagram  for  the  TEDS 
program.  The  FORTRAN  codes  for  the  subroutines  SPIREL  and 
STAN  are  reproduced  in  this  appendix  because  they  are  the 
important  analysis  routines  written  specifically  for  use 
in  this  study.  SPIREL  contains  the  spherical-log  spiral 
curve  fitting  algorithm  and  STAN  produces  the  spiral  analy¬ 
sis  grid  and  calls  the  Fourier  analysis  subroutine  used 
to  quantify  the  spiral  patterns  of  the  cloud  bands. 
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Figure  Al .  Functional  flow  diagram  for  the  TEDS  program. 

The  program  flow  is  down  the  page,  with  horizontal  excur¬ 
sions  where  indicated. 
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COMMON/SCAT/  I  CAT < 256 ) , P H I  X ( 1 00 > . A ( 100) ,B(40> ,C<40) ,0(40) »E(80) 

COMMON /BUFR2/  JOB , I F I LL (  256  ) , I B <  256  ) , NAME  1 (  10) 

COMMON/GF/  I  BUF ( 256 ) 
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COMMON /SV/  I  FLAG , I  HEM , PH  IS ( 20 ) ,XLAMS<  20) .COTALS! 20) 
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EQUIVALENCE  ( AMU ( 1 ) , D ( 1 >  ) 


PI  =  3.1415926535898 
DTR  *  0.0174533 


; DEGREES  TO  RADIANS 


TYPE 

TYPE " WA I T . CACULATIONS  IN  PROGRESS" 

TYPE 


DO  11  1*1,40 
MUO( I  )  *  0.0 
NUO( I >  *  0.0 
AMU (  I  )  *  0.0 
CONTINUE 


XNOPTS  *  F  LOAT ( NOPTS  > 


LONGITUDE  CONVENTION  s  1-360  WITH  WEST  1-180  AND  EAST  181-360 


DO  101  1=1, NOPTS 

LAMBDA! I  )  =AMOD < ( 360 . -LAMB DA( I ) ) , 360 . ) 
IF  (LAMBDA! I ) .EQ.0. >  LAMBDA ( I ) =360 . 
CONTINUE 


XW=AMOD ( ( 360. -XW1 ) , 360. ) 
X E  =  AMOD ( (360.-XE1 ) ,360.  ) 
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DO  12  1*1 , NOP  TS 
LAMBDA ( I >  *LAMBDA ( I >  *DTR 
PHI  ( I  )  *  P  H I ( I >*DTR 
12  CONTINUE 


RMSP  *  1.E70 
SIGN  =  -1 . 

YICOUNT  *  0.0 
XUCOUNT  *  0.0 
COTALFP  *  0.0 
ALEPH  *  0.0 
MU0P  *  0.0 
NU0P  *  0.0 
MUAV  *  0.0 
AVER  *  0.0 
PHIC  *  0.0 
CP H I  *  0.0 
LAMBDAC  =  0.0 
CLAMBDA  *  0.0 

PH  10  *  PHI ( 1 > 

LAMBDA0  =  L AMBDA( 1 > 

DRY  *  (YN  -  YS ) / 6 . 

DRX  =  (XE  -  XW)/6. 

XC  *  XW  +  (DRX  *  3. ) 

YC  *  YS  +  (DRY  *  3. ) 

I  HEM  -  1  ; SET  FLAG  TO  NORTHERN  HEMISPHERE 

IF  ( (PHI0.LT. 0.  ) .AND. (PHI (NOPTS > .LT.0. > )  I  HEM  -  0  ; SOUTHERN  HEMISPHERE 

IF  ( IHEM.EQ.0)  SIGN  =  1 . 


DO  61  I J  *  1,9 


DO  1661  I  *  1,5 
XI  *  FLOAT ( I  ) 

PHIC  *  YC  -  ( 2 . *DR Y )  ♦  ((XI  -  1.)  *  DRY) 
PHIC  =  PHIC  *  DTR 


DO  1661  0  *  1,5 
XO  *  F L OAT ( 0  > 

LAMBDAC  *  XC  -  (2.*DRX>  +  ((XO  -  1.)  *  DRX) 

LAMBDAC  *  LAMBDAC  *  DTR 

CALL  GTLOC  (NU0,  MU0 ,  PHIC,  LAMBDAC,  PHI0,  LAMBDA0 ) 

AMU0  -  ALOG(TAN(MU0*0.5) ) 

SUM  1  ■=  0. 

SUM2  *  0. 

SUM3  »  0. 

SUM4  *  0. 

NUP  =  NU0 
C 
C 

DO  59  K  *  1 ,  NOPTS 

CALL  GTLOC (  NUO(K),  MUO(K),  PHIC,  LAMBDAC,  PHI  (  K ) ,  LAMBDA ( K ) ) 
IF  (MUO(K) .GT.MU0.AND.K.GT. (NOPTS/4) )  GO  TO  1661 
IF  (ABS(NUP-NUO(K) ) .GT.PI )  NUO(K)  -  NUO(K)-(SIGN*2.*PI ) 
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AMU(K)  =  ALOG(TAN(MUO(K)*0.5) ) -AMU0 
SUM1  =  SUM  1  ♦  NUO ( K ) 

SUM2  =  SUM2  ♦  (NUO(K)  *  AMU ( K ) ) 

SUM3  =  SUM3  +  AMU ( K  > 

SUM4  *  SUM4  +  ( AMU ( K ) * AMU ( K ) ) 

NUP  *  NUO(K) 

59  CONTINUE 
C 
C 

CR=(XN0PTS*SUM4  )-(SUM3*SUM3)  ;  DENOMINATOR 

IF  (CR  .EQ.  0. )  COTALF  =  1 . E70 
IF  (CR  .EQ.  0.  )  NU0  =  1 . E  70 
IF  (CR  .EQ.  0. >  GO  TO  567 
COTALF  * ( (SUMI*SUM3 )-<XN0PTS*SUM2) >/CR 
NU0=(  (SUM1*SUM4  )-(SUM2*SUM3)  >/CR 
567  CONTINUE 
C 

RMS  =  0. 

MUSUM  =  0. 


DO  58  K  *  1 , NOP TS 
MUSUM  =  MUSUM  +  MUO(K> 

NUE  =  NU0  -  ( AMU ( K  >  *  COTALF) 

RMS  =  ((NUE  -  NUO  (  K )  >  *  (NUE  -  NUO(K>>)  «■  RMS 
58  CONTINUE 


RMS  =  SQRT(RMS/XNOPTS ) 


IF  (RMS  .GE.  RMSP )  GO  TO  1661 
AVER  =  RMS/SQRT(XNOPTS  > 

MUAV  *  MUSUM/XNOPTS 
RMSP  =  RMS 

COTALFP  *  ABS ( COTAL  F > 

MU0P  =  MU0 
NU0P  *  NU0 
Y I  COUNT  *  FLOAT ( I  ) 

XUCOUNT  -  FLOAT(U) 

CPHI  =  PHIC 
CLAMBDA  =*  LAMBDAC 

1661  CONTINUE 

YC  s  YC  -  (  2 .  ’•'DR  Y )  +  ((YICOUNT  -  1.)  *  DRY) 
XC  =  XC  -  ( 2 .  *DRX )  +  ( ( XUCOUNT  -  1.)  *  DRX ) 
DRY  =  DR Y*0 . 25 
DRX  =  DRX  *0 . 25 

61  CONTINUE 


CPHI  =  CPHI/DTR 
CLAMBDA  =  -CL AMBDA/DTR 

IF  (CLAMBDA.LT. -180. >  CLAMBDA *CLAMBDA  +  360. 

ALEPH  *  ATAN(  1 . /COTALFP  ) 

RADER  *  (MUAV* AVER > / ( COTAL F P * DTR ) 
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c 


TYPE 

TYPE "RMS  LOOP 

TYPE 

TYPE 

TYPE 

TYPE 

TYPE M PH  I C . 

TYPE  H  LAMBDAC .  . 
TYPE 

TYPE"DRX . 

TYP  E " DR Y . 

TYPE 

TYPE " COTAL  F . . . 

TYPE "ALFA . 

TYPE 

TYP  E " N  U0 . 

TYP  E "  MU0 . 

TYPE 

TYPE "RMS . 

TYP  E "AVER . 

TYPE 

TYPE "RADER - 

TYP E " I N  NM - 


COMPLETE;  GOOD  LUCK!  " 


"  ,  C  P  H  I 
M , C  LAMBDA 

”  ,  ( DRX * 4  .  ) 

" , ( DRY*4 .  ) 

" , COTAL  F  P 
" , ( ALEPH/DTR ) 

" , ( NU0P/DTR ) 

" , (MU0P/DTR ) 

" ,  RMS P 

“ , ( AVER /DTR  ) 

" , RADER 
H ,  (  RADER*60 . ) 


WRITE  (18,1000) 

WRITE  (18,2000)  CP  H I , CL AMBDA 
WRITE  (18,3000)  ALEPH/DTR ,COTALFP 
WRITE  (18,5000)  NU0P/DTR ,MU0P/DTR 
WRITE  (18,6000)  RMSP , AVER/DTR 
WRITE  (18,7000)  RADER ,RADER*60. 


TPC  =*  2  .  *  P  I  *  (COS(ALEPH)  *  COS(ALEPH)) 
NUMIN  =  NU0P  +  (SIGN  *  0.05  *  TPC) 

RANGE  *  -SIGN  *  (TPC  +  (4.  *PI>) 

PHIC  *  CP  H I  * DTR 

LAMBDAC  -  AMOD ( ( 3  60 . -CL AMBDA  > ,360.  ) *DTR 
NPRS* 100 

XNPRS  *  FLOAT ( NPRS ) 


CALCULATE  THE  FITTED  SPIRAL 

DO  667  K  *  1  ,  NPRS 
XK  *  FLOAT(K) 

NUEX  *  (  ( XK/XNPRS )  *  RANGE)  +  NUMIN 

MUEX  =  2.*ATAN(TAN(MU0P*0.5)*EXP( (S IGN/COTALF  P )*(NUEX-NU0P ) ) ) 
CALL  LTOGEO  (NUEX,  MUEX,  PHIC,  LAMBDAC,  P H I X ( K ) , LAMB DAX ( K ) ) 
PHIX(K)  -  P  H I  X ( K  )  / DTR 
LAMB  DAX ( K )  *  -L AMBDAX ( K ) /DTR 

IF  (LAMBDAX(K) .LT.-180.  )  LAMB DAX ( K  )  = LAMB DAX ( K )  +  360. 

667  CONTINUE 


C 


DRAW  THE  FITTED  SPIRAL  AFTER  CONVERTING  TO  SCREEN  COORDINATES 
ICOLOR  *  254  ;  GREEN 

CALL  DWSPIR  (PHIX,LAMBDAX, IB, IBUF, NPRS, ICOLOR) 
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c 

c 

c 

c 

c 

c 


c 

c 


c 

c 

c 


c 

c 

c 

c 

999 


C 

C 

C 

1000 
2000 
3000 
5000 
6000 
7000 
1  1000 
C 


DRAW  THE  CENTER  POINT 

I0PT  *  2  ; EARTH  TO  IMAGE  COORDINATES 

CALL  EZNAV  ( I  OPT , I ER , I B , I RAS , I P I X , CP H I , CL AMBDA , 1 ) 

IF  ( IER . NE .0)  GO  TO  999 


ICOLOR  *  256 

CALL  AMOP  (  IPIX-3, IRAS-4) 

CALL  STEXT  ("  +  ",1,1, ICOLOR  *1,1) 
CALL  WLT  ( IBUF, 256,0) 

CALL  XMT  (0) 


TYPE 

ACCEPT"DO  YOU  WISH  TO  SAVE  THIS  SPIRAL’S  DATA  7  ” , IOPT 

IF  ( IOPT . NE . 1 )  GO  TO  999 

I  FLAG  »  I F  LAG  +  1 

P  H IS ( I  FLAG  )  *  PHIC 

XLAMS ( IFLAG)  »  LAMBDAC 

COTALS  < I F  LAG )  =  COTALFP 

WRITE  (18,11000) 


CALL  OVEXIT  < “ OSP IREL " , I  ERR > 

IF  (IERR.NE.l)  TYPE "OSP IREL  EXIT  ERROR 


“ , I  ERR 


FORMAT! ////, 21X , "*«**« 
FORMAT(  IX  ,  ''CENTER  LATITUDE 
FORMAT! IX , "ALFA 
FORMAT! IX , " NU0 
FORMAT! IX, "RMS  ERROR 
FORMAT! IX , “RADIAL  ERROR 
FORMAT! /, 1 IX , "THIS  DATA  WAS 


SPIRAL  OUTPUT  *»««*”,//) 

:  ”  ,E14. 8, 8X, "CENTER  LONGITUDE 
:  "  ,E14.8,8X, "COTAN(ALFA) 

:  " , E 1 4 . 8 , 8X , "MU0 

:  " .E14.8.8X. "AVERAGE  ERROR 
:  " , El  4 . 8 , 8X  ,  "  IN  NM 

SAVED  FOR  THE  SP I RAL -FOUR  I ER 


:  "  ,  E 1  4 . 8  ) 
i  *  ,  E  1  4 . 8  ) 
:  "  ,  E  1  4 . 8  ) 
:  "  ,  E  1  4 . 8  ) 

:  "  ,  E  1  4 . 8  ) 
ANALYSIS . " ) 


END 
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c 

c 


& 


c 

c 


c 

c 

c 

c 


c 

c 


c 


40 


c 

c 

c 


c 

c 


1 

c 

c 


c 

c 

5 


SUBROUTINE  STAN 
OVERLAY  OSTAN 


COMMON /PR/  DATE ,FRMTIMEtIHI , ILO.IRES 
COMMON/HA/  UN ( 100) , NOPTS 

COMMON/SV/  I  FLAG,  I  HEM  ,  P  H  I  S  (  20  >  ,  X  LAMS  (  20 )  tCOTALS(20>  TBAC., 
COMMON /SCAT/  I  CAT (  25  6  ) , PH  I ( 100) , XLAMBDA! 100 ) , I P I X ( 1 00 > , IRAS( 100> , 
F ILL ( 57 ) , AMP ,M90,MAXH , V3  „  „ _ 

COMMON /F/  AMPM,AMP8 , TVAR5 , HEMVM ,V32,TVAR8,V3M, X90M  tTHLD8tSD2 
COMMON /HARM/  COTALF , XMU0 , XNU0 , PH IC , XLAMBDAC 
COMMON /BUFR 1 /  IER,IERR,ILINE(256),IOPT, NTYPE 
COMMON /BUFR2/  JOB,IC(256),IB(256>, NAME  1(10) 

COMMON/GF/  I BUF (  256  ) 


REAL  NU0 , MU0 , LAMBDA! 100) 

REAL  NU , MU 

REAL  SD ( 20 ) , BAR ( 20 ) 

INTEGER  MAX ( 20 ) 

EQUIVALENCE  ( NU0,XNU0) , !MU0,XMU0> , l XLAMBDA! 1 ) , LAMBDA! 1 ) ) 


PI  =  3.1415926535898 
DTR  =  0.0174533 


NOPTS  *  20 
NSECT  =  20 


;  NUMBER  OF  HARMONICS 
; NUMBER  OF  ORTHOGONAL  SLICES 


DO  40  I»1 ,20 
MAX! I  )  =  0 
SD! I )  -  0.0 
BAR! I )  =  0.0 
CONTINUE 


X90M  -  V3M  *  HEMVM  -  AMPM  *  0.0 

XLAMBDAC  =  -XLAMBDAC 

IF  ! XLAMBDAC . LT . - P I >  XLAMBDAC  = 

NOPTS  =  NOPTS*2 

ICOLOR  «  256  ♦ RED 


( 2 . *P I )  ♦  XLAMBDAC 


DO  1  1=1,100 
UN! I )  *  0.0 
CONTINUE 


CALL  CLOSE  !  1  ,  IER > 

CALL  OPEN  !  1  , NAME  1,2, 512, IER) 

IF  !  IER.EQ. 1 )  GO  TO  5 

TYPE  "CANNOT  OPEN  IMAGE  FILE,  IER  *  MER 
IER  =  1 
GO  TO  999 


ALEPH  »  ATANt 1 . /COTALF ) 

TPS  =  2.*PI*!SIN!ALEPH)*SIN!ALEPH)> 
TP  C  =  2.*PI*!COS!ALEPH)*COS!ALEPH) ) 

AMU0  =  TAN ! MU0*0 . 5  > 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

X 

X 

c 

c 

c 


1000 

2000 

5000 


c 

c 

c 

c 


c 

c 


20 


c 

c 

c 


c 

c 


XNOPTS  =  FLOAT ( NOPTS ) 
XNSECT  *  F  LOAT ( NSECT+ 10 ) 
SIGN  =  -1 . 

IF  ( IHEM.EQ.0)  SIGN  *  1 . 


IH  =  I  F  I  X ( FRMT I  ME ) 

WRITE  TO  THE  PARAMETER  FILE: 


NAME  1  *  TROPICAL  CYCLONE  NAME 

DATE  =  YEAR  AND  DATE  OF  TROPICAL  CYCLONE 

IH  *  TIME  OF  IMAGE  OF  TROPICAL  CYCLONE 

I  h  I  =  UPPER  LIMIT  DURING  THRESHOLDING 

Tin  *  LOWER  LIMIT  DURING  THRESHOLDING 

PHIC  =  LATITUDE  OF  SPIRAL  CENTER  USED  IN  ANALYSIS 

LAMBDAC  =  LONGITUDE  OF  SPIRAL  CENTER  USED  IN  ANALYSIS 

XNUf  -  ROTATION  ANGLE  OF  SPIRAL  USED  IN  ANALYSIS 

XMU0  *  SCALE  FACTOR  OF  SPIRAL  USED  IN  ANALYSIS 

ALEPH  *  INFLOW  ANGLE  OF  SPIRAL  USED  IN  ANALYSIS 


F6.0 

12 

13 

13 

F6.2 
F7.2 
F6.2 
F  5 . 2 
F  5 . 2 


WR I TE ( 2  1  , 4000 )  NAME  1 
WRITE (21 ,3000)  DATE, IH, I 

XLAMBDAC  *  AMOD (  <  <2.*PI > 


HI , I LO, PH IC/DTR, XLAMBDAC /DTR, 
-XLAMBDAC) , (2.*PI > ) 


XNU0/DTR, 


XMU0/DTR ,ALEPH/DTR 


WRITE  (18,5000) 

WRITE  (18,1000) 

WRITE  (18,2000)  (0-l,J*l,21) 

FOR MAT ( IX , "WAVENUMBERS " ,/ ) 

FORMAT (21(2X,I2,2X),/) 

FORMAT( ///,30X, "FOURIER  COEF  F I C I  ENTS " , / ) 


DO  10  1*1,  NSECT 

XI  *  FLOAT ( I-l )/XNSECT 


DO  20  0*1  , NOPTS 

XO  *  FLOAT ( 0 ) /XNOPTS 

EX2  *  EXP ( -1 . *TPC*< 1 . +XI )/COTALF ) 

EX 3  *  EXP(COTALF*TPS*XU) 

MU  *  2. *ATAN(AMU0*EX2*EX3) 

NU  =  NU0-(SIGNM  (TPC*( 1  ,+XI  )  )  +  (TPS*XU  )  ) ) 

CALL  LTOGEO  (NU, MU, PHIC, XLAMBDAC, PHI (J)  ,LAMBDA(0)  ) 
PHI (0 ) =  P  H I (0 ) /DTR 
LAMBDA  (  0  )*  -LAMB  DA(  0  )  /DTR 

IF  (  LAMBDA  (  0  ).  LT  1  80  .  )  LAMBDA  (  0  )  *360  .  •♦■LAMBDA(  J  ) 
CONTINUE 


IOPTG  =  T!  THE  DATA  ;  EARTH  TO  IMAGE  COORDINATES 

CALL  EZNAVI IOPT, IER, IB, IRAS, I P IX , PH  I .LAMBDA, NOPTS) 
IF  ( IER.NE.0)  GO  TO  999 


UNBAR  »  0.0 
SMAX  =  0.0 
XSD  =  0.0 
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c 

c 


c 

35 


30 

C 

C 

C 


90 

C 

c 


c 

c 


100 

c 

c 


DO  30  0  =  1 , NOP TS 

CALL  RDBLK( 1 , I R AS ( 0  ) , I  CAT , 1 , I ER  > 

IF  (  IER.EQ. 1 >  GO  TO  35 

TYPE " CANNOT  READ  IMAGE  FILE,  I ER  «  ,IER 
I  ER  *  1 
GO  TO  999 


I  HOLD  *  I  CAT ( (IPIX(J)/2)  +  l) 

IF  ( MOD ( I P I  X ( J ) ,2) .EQ.0)  I  HOLD 
IF  < MOD ( I P IX ( 0 ) , 2 ) . NE .0 )  I  HOLD 
UN ( 0  >  *  FLOAT(IHOLD) 

UNBAR  *  U N { J  )  ♦  UNBAR 

CALL  APOINT  ( I P I  X ( 0 ) , I RAS ( 0 ) , I  COLOR > 
CONTINUE 


ISHFT ( I  HOLD, -8) 
I  HOLD .AND . 377K 


UNBAR  *  UNBAR /XNOPTS 
DO  90  0*1  ,  NOP  TS 

XSD  *  ( (UNBAR-UN( J ) )* ( UNBAR-UN < 0 ) ) )  +  XSD 
CONTINUE 


SMAX  *  UNBAR  -  (  1  .  5 -'SORT <  XS D /XNOPTS  >  ) 
SMAX  *  A  I  NT ( SMAX ) 

UNBAR  *  0.0 
XSD  *  0.0 


DO  100  0=1 , NOPTS 
IF  ( UN ( 0 > . LT . SMAX )  UN(0)  *  SMAX 
UNBAR  *  UN ( 0  >  +  UNBAR 
CONTINUE 


UNBAR  *  UNBAR /XNOPTS 


C 


110 

C 


C 

C 


C 

C 

C 

C 


c 

c 


DO  110  0=1 , NOPTS 

XSD  =  ( ( UNBAR-UN(0 ) )*( UNBAR-UN(0 ) ) )  ♦  XSD 
CONTINUE 

SD  (  I  )  18  SQRT(XSD/XNOPTS) 

BAR ( I >  *  UNBAR 


MAX ( I >  *  I F I X ( SMAX  ) 
IRMAXT  =  I F I  X ( SMAX  > 

IF  ( IRMAXT. GT. 178)  MAXT 
IF  ( IRMAXT. LE. 178)  MAXT 
MAXT  =  MAXT  -  273 


420  -  IRMAXT 
(  662  -  IRMAXT  > /2 


CALL  WLT  ( IBUF, 256,0) 
CALL  XMT  (0) 

CALL  HARMON 


X  90M  =  F  L  OAT ( M90 ) +X  90M 
V3M  =  V3 ♦ V3M 

HEMVM  *  F  LOAT( MAXH )  +  HEMVM 
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ooxxoooooooooo 


AMPM  =  AMP  +AMPM 
IF  ( I .EQ. 16)  AMP 8  =  AMP 
IF  (  I  .  EQ . 4  )  V32  =  V3 
C 

10  CONTINUE 

C 

TVAR5  *  SD( 10)*SD( 10) 
TVAR8  *  SD(  1  6  )  *SD (  16) 
THLD8  *  MAX ( 16) 

SD2  =  SOU) 

X90M  *  X90M/20. 

V3M  *  V3M/20. 

HEMVM  =*  HEMVM/20. 

AMPM  *  AMPM/20. 


WRITE  TO  THE  PARAMETER  FILE: 


BAR (  I  ) 
MAX ( I  ) 
SD  (  I  ) 


*  AVERAGE  VALUE  OF  SLICE 
=  MINIMUM  IRCOUNT  USED  IN  ANALYSIS 
=  STANDARD  DEVIATION  OF  SLICE 


F  8 . 3 
13 

F  5 . 2 


WRITE  (21,7000)  ( BAR ( I ) , I s 1 , 20 ) 

WRITE  (21,6000)  ( MAX (  I  )  ,SD (  I ) , I  =  1 , 20) 


X3000 
X  4000 
X  6000 
X  7000 

C 


F  OR MAT (F6.0,I2,2I3,F6.2,F7.2,F6.2,2F5.2) 
FORMAT? "♦DECK” , IX, 10A2) 

FORMAT? 10 (  I3,F5.2)  ) 

FORMAT  ( 10F8.3,/, 10F8.3) 

CALL  FCST 


999  CALL  CLOSE  (  1  ,  I  ERR  ) 

CALL  OVEXIT  ( " OSTAN " , I  ERR ) 

IF  (IERR.NE.l)  TYPE " OSTAN  EXIT  ERROR  :  ",IERR 
END 
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APPENDIX  B 


GOVERNING  EQUATIONS 

The  slope  of  a  line  on  the  surface  of  a  sphere  is 
given  by 


A  =  w  s^n  0 


(Bl) 


SLOPE 


or 


(B2) 


A 


where  <)>  is  longitude  and  0  is  colatitude.  By  definition 
A  is  constant  for  a  spherical- log  spiral.  Figure  Bl  shows 
the  relationship  between  the  longitudinal  coordinate  <J> 
and  0,  the  colatitude. 

Separating  the  variables  in  (B2)  and  integrating  from 
the  outermost  point  °f  a  spiral  in  towards  the 

csnter ,  we  get  the  equation  for  a  line  of  constant  slope 
on  a  sphere 


4> 

/  d<J>  =  A  / 


and  we  obtain 


<j>  -  <t> 


o 


Jin  tan  (-^-)  1 
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Figure  Bl.  Spherical  polar  coordinates  showing  the  rela¬ 
tionship  between  <{>  and  0 .  r  is  constant  on  the  surface 
of  a  sphere.  4>Q  and  0o  are  the  unit  vectors. 
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or 


A  £n 


tan(|-) 

— % 

tan  (y) 


Renaming  the  variables: 


(B3) 


v  =  <J>  the  azimuthal  coordinate 
and  y  =  6  the  radial  coordinate 


and  defining  the  constant,  A  =  -  cot  a,  where  a  is  the 
spiral  inflow  angle,  (B3)  becomes 


or 


cot  a 


tan(^) 

,  /o, 

tan  (y) 


v 


and 


v 


o 


cot  a  £n 


tan  (^-) 

rx 

tan  (y) 


(B4) 


=  2  arctan  (tan(y)  exp[-  tan  a(v-VQ)]} 


(B5) 


Equation  (B5)  is  the  equation  for  a  logarithmic  spiral  de¬ 
creasing  with  increasing  azimuth  and  having  a  constant  in¬ 
flow  angle  (a)  on  the  surface  of  a  sphere.  Equation  (B4) 
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is  in  the  form  of  an  equation  for  a  line: 


y  =  mx  +  b 


(B6) 


where : 

y  =  v, 

tan  ( j) 

x  =  £n  -  , 

yo 

[tan  (-2-)  J 

b  =  v  ,  and 
o 

m  =  -  cot  a. 

Therefore,  all  points  along  a  spiral,  y^  and  v^,  can  be 
plotted  on  a  cartesian  graph  with  abscissa  £n [ (tan (y/2) ) / 
(tan(yQ/2) ) ]  and  ordinate  v;  and  a  spherical-log  spiral 
would  be  represented  by  a  straight  line  with  slope  -  cot  a. 
To  map  any  point  on  the  surface  of  a  sphere  into  the  spiral 
coordinates  the  latitude-longitude  representation  of  the 
point  must  be  converted  to  y  and  v,  the  great  circle  dis¬ 
tance  between  the  point  and  the  spiral  center  and  their 
angular  separation,  respectively. 

In  order  to  specify  the  spiral  coordinates  the  spiral 
center  is  needed  and  is  furnished  by  convergent  search 
routine  which  returns  a  best  fit  center  in  a  least  squares 
sense.  The  problem  of  fitting  the  best  least  squares  line 
to  a  collection  of  data  involves  minimizing 
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N  2 
l  [Y.  -  (mX.  +  b) ] 

i=l  1 


with  respect  to  m  and  b. 

Equation  (B6)  yields  the  normal  equations,  which  must 
be  satisfied  for  a  minimum  to  occur. 

N  ? 

m  I  Xz  +  b 

•  T  1 

1  =  1 


and 


N 

m  7  X .  +  bN 
i=l 


where : 


m  =  -  cot  a 


X. 

l 


in 


r  /i,  1 

tan  (~2~) 
Ltan  (-^-)  J 


The  solution  to  (B7)  is  the  following: 
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m  = 


N 


N  N  N 

l  X  Y  -  I  x  l  Y± 

i=l  1  1  i=l  1  i=l 

N  ,  N  2 

N  l  X*  -  (I  X  ) 

1=1  1=1 


and 


(B8) 


b  = 


N 

l 

i=l 


X  , 


N 

l 

i=l 


Y.  - 


N  N 

I  X  Y  l  X 
i=l  1  1  i=l  1 


N 


N  2 

Ji  <  - 


N 


(  l  x.) 

i=l  1 


From  the  estimated  center  and  yQ/  vq  and  a  are  calcu¬ 
lated  using  the  above  linear  least  squares  fit.  Using  the 
resulting  spiral  parameters  the  RMS  error  is  calculated 
for  the  input  points  and  used  as  the  criterion  for  con¬ 
verging  on  a  more  accurate  guess  of  the  spiral  center . 

The  convergent  search  routine  divides  the  area  esti¬ 
mated  to  contain  the  spiral  center  into  a  six-by-six  grid. 
Each  of  the  twenty-five  interior  grid  points  are  used  as 
sample  spiral  centers  and  the  RMS  error  calculated.  The 
grid  point  with  the  smallest  RMS  error  is  then  focused  on 
as  the  new  grid  center  and  the  grid  size  divided  by  six. 
The  process  is  iterated  nine  times,  with  the  final  center 
the  best  fit  (smallest  RMS  error) . 

The  equation  for  an  orthogonal  spiral  is  just  (B6) 
with  m  replaced  by  - 
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y 


=  -  -  x  +  b 

m 


or,  in  spiral  notation 


tan  (~) 
tan  (“2“) 


o 


v 


cot  a 


+  v 


(B9) 


APPENDIX  C 


DESCRIPTION  OF  CYCLONES  OF  DEPENDENT  DATA  SET 

1 .  DOLORES 

Hurricane  Dolores  was  the  sixth  tropical  cyclone  of  the 
Eastern  North  Pacific  (ENP)  1979  season.  The  initial  dis¬ 
turbance  was  spawned  west  of  the  Central  American  coast  on 
14  July  over  29 °C  water.  The  disturbance  was  upgraded  to 
a  tropical  depression  by  0600  GMT  17  July  and  further  to 
a  tropical  storm  by  1800  GMT  17  July  as  the  winds  increased 
above  35  kt.  By  the  next  day  the  winds  had  increased  to 
above  65  kt  and  Hurricane  Dolores  was  born.  At  this  time 
it  changed  course  from  westerly  to  north-westerly,  and 
began  a  rapid  intensification.  An  eye  became  visible  on 
satellite  photographs  by  1800  GMT  19  July  and  by  0000  GMT 
21  July  Dolores  had  reached  her  maximum  intensity  with  winds 
of  over  105  kt.  Accelerating  around  an  upper  level  high- 
pressure  area  centered  over  Mexico,  Dolores  moved  to  a  posi¬ 
tion  thirteen  degrees  west  of  the  Baja  California  coast  and 
began  to  weaken  over  cool  water.  By  1800  GMT  22  July, 
Dolores  was  downgraded  to  a  tropical  depression  and  finally 
dissipated  with  a  cloud  mass  lingering  in  the  Pacific  for 
three  days . 

2.  ENRIQUE 

Hurricane  Enrique,  the  seventh  tropical  cyclone  of  the 
ENP  1979  season,  developed  as  a  tropical  disturbance  near 
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10°  N,  108°  W  on  17  August.  At  1800  GMT  that  day  the  dis¬ 
turbance  was  upgraded  to  a  tropical  depression  and  by  0000 
GMT  18  August  was  upgraded  again  to  Tropical  Storm  Enrique. 

At  0000  GMT  19  August  the  west-northwestward  moving  Enrique 
was  declared  Hurricane  Enrique  and  began  to  show  an  eye  on 
satellite  photographs.  The  maximum  winds,  70  kt,  decreased 
to  55  kt  by  0600  GMT  21  August  as  then  Tropical  Storm  Enrique 
moved  over  low  clouds  and  cool  water .  However,  by  1700 
GMT  that  same  day  the  storm  was  upgraded  to  hurricane  status 
again  as  the  winds  increased  above  70  kt.  Enrique's  devel¬ 
opment  was  now  rapid  as  it  reached  its  maximum  intensity  at 
1200  GMT  22  August  with  winds  of  125  kt.  Enrique  began  to 
weaken  as  it  moved  over  cool  water  and  low  clouds,  and  by 
1200  GMT  23  August  was  downgraded  back  to  tropical  storm 
status.  By  0000  GMT  24  August  the  storm  was  further  down¬ 
graded  to  a  tropical  depression  and  weakened  rapidly,  dissi¬ 
pating  over  the  next  two  days. 

3 .  DAVID 

David  was  the  fourth  hurricane  in  the  Atlantic  1979 
season.  It  organized  in  the  mid-tropical  Atlantic  and  was 
declared  a  tropical  depression  by  0000  GMT  25  August. 
Tracking  westward  around  the  southern  portion  of  the  Ber¬ 
muda  High,  the  disturbance  was  upgraded  to  a  tropical  storm 
by  1800  GMT  25  August.  The  storm  continued  its  steady  in¬ 
crease  in  intensity  and  was  upgraded  to  hurricane  strength 
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by  1800  GMT  26  August.  Hurricane  David  continued  its  de¬ 
velopment  and  west-northwestward  movement  and  reached  the 
eastern  coast  of  Hispaniola  on  31  August  with  a  central 
pressure  of  approximately  925  mb  and  estimated  winds  in 
excess  of  150  kt.  Flooding  caused  by  the  excessive  rain¬ 
fall  from  David  was  responsible  for  over  2000  deaths  and 
over  1  billion  dollars  damage  on  Hispaniola  (Hebert,  1980)  . 
David  continued  its  clockwise  track  and  entered  Florida  by 
early  3  September  with  winds  in  excess  of  80  kt.  David 
weakened  slowly  over  land,  causing  15  deaths,  damages  of 
320  million  dollars  and  evacuations  of  400,000  people 
(Hebert,  1980) .  David  was  finally  downgraded  to  tropical 
storm  status  over  South  Carolina  at  1800  GMT  4  September. 
Tropical  storm  David  was  further  downgraded  to  depression 
stage  on  6  September  as  its  winds  fell  below  35  kt  and  cen¬ 
tral  pressure  increased  to  992  mb.  Strong  winds  from  the 
remains  of  Hurricane  David  continued  to  persist  along  the 
Atlantic  seaboard  all  the  way  into  the  Canadian  provinces. 

4 .  FREDERIC 

Hurricane  Frederic  started  life  as  an  African  tropical 
wave  on  27  August  and  turned  out  to  be  the  sixth  hurricane 
of  the  1979  Atlantic  season.  The  wave  continued  to  develop 
and  reached  tropical  depression  stage  by  1800  GMT  28  August; 
the  disturbance  became  Tropical  Storm  Frederic  by  0000  GMT 
30  August  while  continuing  on  its  western  track.  At  1800 
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GMT  31  August  Tropical  Storm  Frederic  became  Hurricane 
Frederic  in  the  same  area  of  the  Atlantic  that  David  had 
done  four  days  before;  an  eye  became  visible  on  the  satel¬ 
lite  imagery.  However,  Frederic  began  to  interact  with 
the,  by  then,  immense  outflow  from  David  and  weakened  back 
to  tropical  storm  status.  Frederic  passed  over  Puerto  Rico 
and  further  weakened  to  depression  stage  by  0600  GMT  6 
September.  The  disturbance  then  turned  northwestward  around 
the  high  and  strengthened  back  to  tropical  storm  intensity 
by  0000  GMT  9  September  off  of  the  southern  coast  of  Cuba. 
Freed  from  the  influence  of  David,  Frederic  regained  hurri¬ 
cane  force  over  western  Cuba  with  a  central  pressure  990  mb 
at  1200  GMT  10  September.  Hurricane  Frederic  continued  to 
develop  and  passed  over  the  Gulf  Coast  of  the  United  States 
with  a  central  pressure  of  946  mb,  and  winds  of  approxi¬ 
mately  115  kt,  at  0300  GMT  13  September.  Frederic  caused 
over  2.3  billion  dollars  damage  in  the  United  States,  mak¬ 
ing  it  the  costliest  hurricane  in  history  for  the  United 
States  (Hebert,  1980) .  Hurricane  Frederic  quickly  filled 
and  was  degraded  to  tropical  storm  status  over  Mississippi 
by  1200  GMT  13  September.  Frederic's  remnants  continued  to 
track  over  the  Midwest,  dissipating  to  depression  status 
by  1200  GMT  14  September  and  exiting  into  Canada  by  0000 
GMT  15  September. 
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5 .  GLORIA 


Gloria  was  the  seventh  hurricane  of  the  Atlantic  1979 
season  and  was  the  first  of  the  season  not  to  affect  a 
landfall.  The  disturbance  started  as  a  well  organized 
African  tropical  wave  and  was  classified  as  a  tropical 
depression  at  1200  GMT  4  September.  The  depression  moved 
into  the  mid-Atlantic  and  became  Tropical  Storm  Gloria  by 
1200  GMT  6  September.  By  0000  GMT  7  September  Gloria  was 
upgraded  to  hurricane  strength  and  started  to  move  north¬ 
westward  in  the  mid-Atlantic.  Gloria  intensified  for  a 
day,  weakened  briefly,  then  regained  intensity  until  it 
reached  its  lowest  central  pressure  of  975  mb  on  1800  GMT 
12  September.  Gloria  then  began  to  fill  and  move  northeast¬ 
ward  and  was  downgraded  to  a  tropical  storm  by  0600  GMT 
15  September  and  eventually  merged  with  an  extratropical 
low  pressure  system  northwest  of  the  Azores  by  late  that 
day . 

6.  IGNACIO 

Hurricane  Ignacio  was  the  twelfth,  of  thirteen,  tropi¬ 
cal  cyclone  of  the  ENP  1979  season.  It  began  as  a  tropi¬ 
cal  disturbance  off  of  the  Central  American  coast,  reached 
tropical  depression  status  by  1800  GMT  23  October  positioned 
near  12°  N,  95°  W.  The  depression  increased  in  intensity 
while  over  warm  water  and  was  declared  Tropical  Storm 
Ignacio  at  1200  GMT  24  October.  The  storm  moved  west  for 
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a  day,  then  northwest,  paralleling  the  Mexican  coast,  and 
was  upgraded  to  Hurricane  Ignacio  by  1800  GMT  26  October  as 
its  winds  increased  to  70  kt.  Ignacio  reached  its  maximum 
intensity  by  2200  GMT  27  October  near  17°  N,  106°  W. 

Ignacio  continued  to  move  northwest  for  a  day  then  recurved 
and  began  to  weaken.  The  hurricane  was  downgraded  to  tropi¬ 
cal  storm  status  by  1800  GMT  29  October  as  the  winds 
dwindled  past  55  kt.  By  1200  GMT  30  October  the  storm  was 
downgraded  to  a  tropical  depression  and  brought  a  landfall 
over  the  Central  Mexican  coast  slightly  north  of  Acapulco. 
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APO  SAN  FRANCISCO  96519 


COMMANDING  OFFICER 
NORDA 

NSTL  STATION,  MS  39529 


COMMANDING  OFFICER  (2) 

NAVWESTOCEANCEN 

BOX  113 

PEARL  HARBOR,  HI  96860 


COMMANDING  OFFICER 
U.S.  NAVOCEANCOMCEN 
BOX  12,  COMNAVMARIANAS 
FPO  SAN  FRANCISCO  96630 


SUPERINTENDENT 
LIBRARY  REPORTS 
U.S.  NAVAL  ACADEMY 
ANNAPOLIS,  MD  21402 
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NAVAL  POSTGRADUATE  SCHOOL 
METEOROLOGY  DEPT. 
MONTEREY,  CA  93943 


COMMANDER 

NAVAIRSYSCOM  (AIR-330) 
WASHINGTON,  DC  20361 


USAFETAC/TS 
SCOn  AFB,  IL  62225 


DFT  R  iww /rc 

APO  SAN  FRANCISCO  96274 


NOAA-NESDIS  LIAISON 
ATTN:  DR.  M.  HELFERT,  CODE  SC2 
NASA-JOHNSON  SPACE  CENTER 
HOUSTON,  TX  77058 


CHIEF,  SCIENTIFIC  SERVICES 
NWS/NOAA,  SOUTHERN  REGION 
ROOM  10E09 
819  TAYLOR  ST. 

FT.  WORTH,  TX  76102 


METEOROLOGIST  IN  CHARGE 
WEA.  SERV.  FCST.  OFFICE,  NOAA 
660  PRICE  AVE. 

REDWOOD  CITY,  CA  94063 


CHIEF 

MESOSCALE  APPLICATIONS  BRANCH 
NATIONAL  EARTH  SAT.  SERV. 

1225  W.  DAYTON 
MADISON,  WI  53562 


LIBRARY 

NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CA  93943 


COMMANDER 

NAVAIRSYSCOM 

MET.  SYS.  DIV.  (AIR-553) 

WASHINGTON,  DC  20360 


AFGL/LY 

HANSCOM  AFB,  MA  01731 


DIRECTOR  (12) 

DEFENSE  TECH.  INFORMATION 
CENTER,  CAMERON  STATION 
ALEXANDRIA,  VA  22314 


DIRECTOR 

NATIONAL  EARTH  SAT.  SERV/SEL 
FB-4,  S321B 
SUITLAND,  MD  20233 


CHIEF,  SCIENTIFIC  SERVICES 
NWS/NOAA,  WESTERN  REGION 
P.O.  BOX  11188, 

FEDERAL  BLDG. 

SALT  LAKE  CITY,  UT  84111 


METEOROLOGIST  IN  CHARGE 
WEA.  SERV.  FCST.  OFFICE,  NOAA 
HONOLULU  INTL.  AIRPORT 
HONOLULU,  HI  96819 


DIRECTOR  (AOML) 

NATIONAL  HURRICANE  RSCH.  LAB. 
1320  S.  DIXIE  HWY. 

CORAL  GABLES,  FL  33146 


COMMANDER  (2) 

NAVAIRSYSCOM 
ATTN:  LIBRARY  (AIR-7226) 
WASHINGTON,  DC  20361 


COMMANDER 

AWS/DN 

SCOTT  AFB,  IL  62225 


DET  4  HQ  AWS/CC 

APO  SAN  FRANCISCO  96334 


ENVIRONMENTAL  SERVICES  DIV. 
OFFICE  OF  THE  JOINT  CHIEFS 
OF  STAFF 
THE  PENTAGON 
WASHINGTON,  DC  20301 


DIRECTOR 

NATIONAL  HURRICANE  CENTER 
NOAA,  GABLES  ONE  TOWER 
1320  S.  DIXIE  HWY 
CORAL  GABLES,  FL  33146 


CHIEF,  SCIENTIFIC  SERVICES 
NWS,  PACIFIC  REGION 
P.O.  BOX  50027 
HONOLULU,  HI  96850 


DIRECTOR,  CENTRAL  PACIFIC 
HURRICANE  CENTER  NWS,  NOAA 
HONOLULU,  HI  96819 


DIRECTOR 

TECHNIQUES  DEVELOPMENT  LAB 
GRAMAX  BLDG. 

8060  13TH  ST. 

SILVER  SPRING,  MD  20910 
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LABORATORY  FOR  ATMOS.  SCI. 

NASA  GODDARD  SPACE  FLIGHT  CEN. 
GREENBELT,  MD  20771 


COLORADO  STATE  UNIVERSITY 
ATMOSPHERIC  SCIENCES  DEPT. 
ATTN:  LIBRARIAN 
FT.  COLLINS,  CO  80523 


CHAIRMAN,  METEOROLOGY  DEPT. 
PENNSYLVANIA  STATE  UNIV. 

503  DEIKE  BLDG. 

UNIVERSITY  PARK,  PA  16802 


UNIVERSITY  OF  MARYLAND 
METEOROLOGY  DEPT. 
COLLEGE  PARK,  MD  20742 


AMERICAN  METEORO.  SOCIETY 
METEOR.  &  GEOASTRO.  ABSTRACTS 
P.O.  BOX  1736 
WASHINGTON,  DC  20013 


BUREAU  OF  METEOROLOGY 
BOX  1289K,  GPO 
MELBOURNE,  VIC,  3001 
AUSTRALIA 


DIRECTOR,  TYPHOON  MODERATION 
RSCH  &  DEV.  OFFICE,  PAGASA 
MINISTRY  OF  NATIONAL  DEFENSE 
1424  QUEZON  AVE. 

QUEZON  CITY,  PHILIPPINES 


COLORADO  STATE  UNIVERSITY 
ATMOSPHERIC  SCIENCES  DEPT. 
ATTN:  DR.  WILLIAM  GRAY 
FORT  COLLINS,  CO  80523 


CHAIRMAN,  METEOROLOGY  DEPT. 
METEORO.  &  SPACE  SCI.  BLDG. 
1225  W.  DAYTON  ST. 

MADISON,  WI  53706 


FLORIDA  STATE  UNIVERSITY 
ENVIRONMENTAL  SCIENCES  DEPT. 
TALLAHASSEE,  FL  32306 


R.  CECIL  GENTRY 
DEPT.  OF  PHYSICS  &  ASTRONOMY 
CLEMSON  UNIVERSITY 
CLEMSON,  SC  29631 


MR.  W.  G.  SCHRAMM/WWW 
WORLD  METEOROLOGICAL 
ORGANIZATION 

CASE  POSTALE  #5,  CH-1211 
GENEVA,  SWITZERLAND 


METEOROLOGICAL  INSTITUTE 
KYOTO  UNIVERSITY 
ATTN:  DR.  R.  YAMAMOTO 
SAKYO,  KYOTO  606,  JAPAN 


CHAIRMAN,  METEOROLOGY  DEPT. 
CALIFORNIA  STATE  UNIVERSITY 
SAN  JOSE,  CA  95192 


DIRECTOR,  REMOTE  SENSING  LAB 
UNIVERSITY  OF  MIAMI 
P.O.  BOX  248003 
CORAL  GABLES,  FL  33124 


UNIVERSITY  OF  HAWAII 
METEOROLOGY  DEPT. 
2525  CORREA  ROAD 
HONOLULU,  HI  96822 


THE  EXECUTIVE  DIRECTOR 
AMERICAN  METEORO.  SOCIETY 
45  BEACON  ST. 

BOSTON,  MA  02108 


DIRECTOR,  JTWC  (3) 

BOX  17 

FPO  SAN  FRANCISCO  96630 


LIBRARIAN,  NEW  ZEALAND 
METEOROLOGICAL  SERVICE 
P.O.  BOX  722 
WELLINGTON,  NEW  ZEALAND 
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